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A quasiparticle pattern advanced in Landau's first article on Fermi liquid theory is adapted to 
elucidate the properties of a class of strongly correlated Fermi systems characterized by a Lifshitz 
phase diagram featuring a quantum critical point (QCP) where the density of states diverges. The 
necessary condition for stability of the Landau Fermi Liquid state is shown to break down in such 
systems, triggering a cascade of topological phase transitions that lead, without symmetry violation, 
to states with multi-connected Fermi surfaces. The end point of this evolution is found to be an 
exceptional state whose spectrum of single-particle excitations exhibits a completely flat portion at 
zero temperature. Analysis of the evolution of the temperature dependence of the single-particle 
spectrum yields results that provide a natural explanation of classical behavior of this class of Fermi 
systems in the QCP region. 



I. PREAMBLE 

Once upon a time at a traditional meeting of nestlings 
of the Migdal school, A. B. began his speech with words 
that stunned the audience. 

Until recent times I was proud that I have never pub- 
lished an incorrect article. But is it actually a true cause 
for pride? Can an experienced mountaineer take pride 
in not having broken any ribs, or a professional motorcy- 
clist, that his legs are still intact? The correct answer is 
no! It just means that these people may have run to the 
best of their abilities, but certainly not more. 

This memorial article, dedicated to a great physicist 
and a great man, is devoted to a problem first discussed 
around 20 years agOii~— The cited works considered the 
possibility of a breakdown of the conventional Landau- 
Migdal quasiparticle patterniri of phenomena observed 
in Fermi liquids (FL), associated specifically with rear- 
rangement of the T = Landau quasiparticle momentum 
distribution np(p) = 0(pf — p)- 

Over the last decade, experimental studies of non- 
Fermi liquid (NFL) behavior of strongly correlated sys- 
tems have extended the frontiers of low-temperature con- 
densed matter physics^— During the same period, a 
number of theorists have engaged in efforts to extend the 
frontiers of FL theory with the aim of explaining this 
anomalous behavior 

There is a famous Migdal Correspondence Criterion for- 
judging new theories, which boils down to this: 

First an foremost, the proposed theory must be able 
to match results of previous, well-tested theoretical de- 
scriptions. It is only of secondary importance that it 
matches the relevant experimental data. Why so? Be- 
cause as a rule, the experiments of burning interest were 
performed yesterday, and therefore the results obtained 
may be flawed, whereas theoretical physics stands as body 
of knowledge created and honed by a multitude over three 



hundred years. 

The extended quasiparticle picture to be reviewed and 
discussed here meets the Migdal Criterion. It reduces 
naturally to the the standard FL picture when deal- 
ing with conventional Fermi liquids. At the same time, 
convincing experimental evidence supporting this exten- 
sion of the Landau-Migdal vision, while present, remains 
scarce — possibly because of the short time the relevant 
experimental programs have been in operationa l 14 ' 29 i 30 
Even so, one cause for the lack of unambiguous corre- 
lations between available experimental data and corre- 
sponding theoretical results may be conceptual (or tech- 
nical) errors made by the theorists in developing the new 
theory, errors that are yet to be exposed. In this case, 
we should not be too upset, since according to Migdal's 
provocative challenge to his assembled group, such errors 
tell us that the struggling theorists, being professionals, 
did, after all, run beyond the breaking point of their abil- 
ities. 



II. ROUTES TO BREAKDOWN OF STANDARD 
FERMI-LIQUID THEORY 

Any theory has own limits of applicability, and FL the- 
ory is no exception. Conventionally these limits are im- 
puted to violation of Pomeranchuk stability conditions. 31 
Such violation is associated with second-order phase 
transitions that are accompanied with jumps of the spe- 
cific heat C(T) and cusps of the spin susceptibility %(T). 
Upon approach to the transition point, spontaneous cre- 
ation and enhancement of fluctuations suppresses the 
value of the z factor that determines the quasiparticle 
weight in the single-particle state at the Fermi surface. 
At the transition point, z vanishes, signaling a breakdown 
of the FL quasiparticle picture ! 32 ' 33 

A different domain where predictions of FL theory 
prove to be fallacious has been discovered and explored 
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during the last decade. This is the regime of the so-called 
quantum critical point (QCP) revealed in experimental 
studies of 2D liquid 3 He and heavy-fermion metals. In 
this domain, the specific heat C(T) and spin suscepti- 
bility x(T), both proportional to the density of states 
N(T) (and hence to the effective mass M*), are found to 
diverge as the temperature T goes to zero4~— 

The presence of a QCP is the hallmark that distin- 
guishes strongly correlated Fermi systems from those 
with weak or moderate correlations. In this respect, 
three-dimensional (3D) liquid 3 He, for which M* remains 
finite at any density, belongs to the class of systems with 
moderate correlations. Its 2D counterpart, on the other 
hand, is assigned to the class of strongly correlated Fermi 
systems based on evidence for the existence of a QCP 
provided by experimental studies of dense 3 He films£~— 

It is well to emphasize that the divergence of the 
specific heat C(T) that is observed in the QCP region 
has little in common with the jumps in C(T) inher- 
ent in second-order phase transitions. Notwithstand- 
ing this fact, theoretical explanations of the failure of 
FL theory in the QCP region have commonly linked its 
apparent breakdown with attendant second-order phase 
transitions! 13 i 32 ~ — However, as they occur in the QCP re- 
gion, such transitions are found to be quite atypical. In 
most cases, even the structure of the order parameters 
remains unknown, and the properties of states beyond 
the QCP defy explanation within the standard scaling 
theory of second-order phase transitions. One is led to 
conclude that the viability of such transitions as triggers 
of the observed NFL behavior is problematic. 

In this article, we pursue another explanation of NFL 
behavior in the QCP region, attributing the breakdown 
of standard FL theory to violation of the necessary sta- 
bility condition for the Landau stated The salient fea- 
ture of the proposed scenario is rearrangement of single- 
particle degrees of freedom driven by topological phase 
transitions (TPT's). This mechanism stands in direct 
contrast to the violation of sufficient Pomeranchuk sta- 
bility conditions, 31 which entail rearrangement of collec- 
tive degrees of freedom in second-order phase transitions. 

Since the topological scenario docs not implicate 
the violation of any Pomeranchuk stability condition, 
the catastrophic suppression of the quasiparticle weight 
posited in the collective scenario for the QCP does not 
take place, and failure of standard FL theory in the QCP 
region must have another explanation. Indeed, the fun- 
damental FL formula 



in a topological rearrangement of the Fermi surface that 
necessitates modification of the standard FL formalism. 



III. REPRISE OF THE STANDARD FL 
QUASIPARTICLE PICTURE 

To set the stage for discussion of the required modi- 
fications, we recall that the heart of the Landau-Migdal 
quasiparticle picture is the postulate that there exists a 
one-to-one correspondence between the totality of real, 
decaying single-particle excitations of the actual Fermi 
liquid and a system of immortal interacting quasiparti- 
cles. There are two facets of this correspondence. First, 
the number of quasiparticles is equal to the given number 
N of particles, a condition expressed as 



r n 

Try n{p)dv = - 



(2) 



where n(p) is the quasiparticle momentum distribution, p 
is the density, Tr implies summation over spin and isospin 
variables, and dv = dp/(2ir) D is a volume element in a 
momentum space of dimension D. Second, all thermody- 
namic quantities, notably the ground state energy E and 
entropy S, are treated as functionals of the quasiparticle 
momentum distribution ?i(p). In particular, the entropy 
S of the real system is given by the ideal- Fermi-gas com- 
binatorial expression 
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V 



-Tr J [n(p) In n(p) + (1 - n(p)) ln(l - n(p))]dv. 

(3) 

In homogeneous matter, a standard variational procedure 
based on Eq. ^ and involving restrictions that impose 
conservation of the particle number and energy, leads to 
the result 4 



i(p) = 1 



(4) 



While this relation resembles the corresponding Fermi- 
Dirac formula for the ideal Fermi gas, the quasiparti- 
cle energy e(p) measured from the chemical potential 
does not coincide with the bare single-particle energy 
e p — P 2 1 (2M) — fi. Instead, it is given by the variational 
derivative 



e(p) = Sn/6n(p) 



(5) 



i{p) = v F (p-p F ) = p F (p - p F )/M* 



(1) 



for the single-particle spectrum e(p) measured from the 
Fermi surface, becomes powerless when the effective mass 
M* diverges. If such a divergence is present, additional 
terms of the Taylor expansion of e(p) must be added to 
the right side of Eq. ([I]) (see Ref. |22| ). giving rise to a 
change of sign of the Fermi velocity. This in its turn sig- 
nals a violation of the necessary stability condition for the 
conventional Landau state (as seen in Sec. V.B), resulting 



of the thermodynamic potential Q — E — /jlN . 

It is significant that the most seminal article cited by 
Landau in his pioneering works^ was authored by A. B. 
Migdal, 36 who had advanced an idea that allows unam- 
biguous determination of the Fermi momentum pp of the 
quasiparticle system. In his numerous papers and books, 
A. B. also adapted the Landau quasiparticle concept to 
finite systems where the phenomenological spectrum (fTJ) 
becomes more complicated due to its dependence on the 
shape of the self-consistent field acting on a quasiparticle 
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in an inhomogeneous medium. Indeed, his famous book 6 
completed the creation of the FL formulation through 
the exposition of a quantitative theory of the behavior 
of Fermi liquids in external fields. This book has be- 
come a touchstone for modern theoretical description of 
the properties of atomic nuclei. Moreover, together with 
A. I. Larkin he developed a quantitative theory of super- 
fluid Fermi liquids^ It is in recognition of these splendid 
achievements that the FL quasiparticle theory is often re- 
ferred to as the Landau-Migdal quasiparticle pattern. 

In principle, the quasiparticle spectrum e(p) can be 
evaluated by means of the relatio n 4 ' 5 ' 38 i 39 



v(p) 



de{p) 
dp 



M 



/(p,Pi 



dn(pi 
dpi 



■dvi, (6) 



stemming from the assertion that the single-particle en- 
ergy e(p) is a functional of the quasiparticle momentum 
distribution n(p). Alternatively, this equation may be de- 
rived from the Galilean invariance of the system Hamil- 
tonian and gauge invariance ! 40 ' 41 From this coincidence 
we infer that the spectrum e(p) is, indeed, a functional 
of n, thereby reinforcing the validity of the basic Lan- 
dau postulate that the ground state energy E and other 
thermodynamic quantities can be treated as functionals 
of n. 

Eq.Q provides a nonlinear integral equation for self- 
consistent determination of the quasiparticle spectrum 
e(p) and the quasiparticle momentum distribution 
treating the Landau interaction function /(p, pi) as phe- 
nomenological input. It is worth noting that the correct 
relation between the function f(p, pi) and the scattering 
amplitude T was first established by Migdal. 42 This func- 
tion turns out to be the so-called w-limit of T(p, p 1; k, w), 
where the energy transfer u) and momentum transfer k 
both tend to zero, but in such a way that k/uj —> 0. 
In this limit, only regular Feynman diagrams contribute 
to / in full force; consequently the interaction function 
/ cannot be evaluated within the FL formalism itself. 
Accordingly, / is commonly specified by a set of phc- 
nomenological parameters, namely the harmonics of its 
Legendre polynomial expansion. 

Following Landau (cf. Eq. (4) in Ref. 0), one postu- 
lates that at T = 0, solutions of Eq. © arrange them- 
selves in such a way that the Fermi velocity vp always 
keeps a positive value. Necessarily, then, the quasipar- 
ticle momentum distribution n(p, T = 0) coincides with 
the Fermi step np{p) — 9{pf~p) appropriate to a nonin- 
teracting Fermi gas. In turn, the relation ([2]) reduces to 
the famous Landau-Luttinger theorem, expressed simply 



P : 



3tt 2 



(7) 



for the case of 3D homogeneous matter. The resulting 
quasiparticle picture, in which a system of fermions is 
treated as a "gas of interacting quasiparticles" 6 works 
flawlessly for conventional Fermi liquids, including 3D 



liquid 3 He, the matter inside atomic nuclei, and the elec- 
tron liquid in alkali metals, where the correlations are 
moderate or weak. However, as already indicated, this 
picture begins to fail upon entry into the QCP region, 
where the effective mass M* diverges. 



IV. THEORETICAL ALTERNATIVES FOR 
DIVERGENCE OF THE EFFECTIVE MASS 

Discrepancies between FL predictions and experimen- 
tal data were first observed in films of liquid 3 He at rel- 
atively high areal density. In stark contrast to the ordi- 
nary FL behavior exhibited by the bulk liquid, the 2D 
data for the spin susceptibility x(T") and Sommerfeld ra- 
tio 7 = C{T)/T of specific heat C(T) to temperature T 
soar upward as T — > 0,— ~— whereas both quantities must 
become constant in FL theory. Analogous behavior has 
been discovered and documented in strongly correlated 
electron systems of solids 1 ^— and in the 2D electron gas 
as realized in MOSFETs.— ~— 

Such a behavior is ascribed to divergence of the effec- 
tive mass M*, first revealed theoretically in microscopic 
calculations of the single-particle spectrum e(p) of the di- 
lute 3D homogeneous electron ga s 48 ' 49 (for the 2D case, 
see Ref. Two different sources for divergence of M* 
in nonsuperfluid Fermi systems can be identified in terms 
of the textbook formula 



M 
~M* 



de° p 



(8) 



where the quasiparticle weight z in the single-particle 
state is given by z = [1 — (<9E(p, s)/de) Q ]~ 1 , with E rep- 
resenting the mass operator. Here and henceforth, the 
subscript indicates that a function of p and e is evalu- 
ated at the Fermi surface. 

As seen from this formula, one way in which a di- 
vergence of M* could arise (first discussed by Doniach 
and Engelsberg 3 -^) is through divergence of the deriva- 
tive (<9E(p, e)/de) at a critical density p c where some 
second-order phase transition occurs. Another option in- 
vokes a change of sign of the sum 1 + (<9E (j?,e)/ 9e°) at a 
critical density . Correspondingly, there are two differ- 
ent scenarios for the QCP. In the first scenario — the most 
popular until recently — it is the energy dependence of the 
self-energy E(p, e) that plays the decisive role. In this col- 
lective scenario, the posited divergence of (dT,(p,e)/de) 
at the points of second-order phase transitions results 
in a vanishing renormalization factor z, and the QCP is 
identified with the end point of the line T/v(p) traced in 
the phase diagram by a second-order phase transition. 
The essence of this scenario is captured in the maximi 3 - 
that in the vicinity of any second-order phase transition, 
"Quasiparticles get heavy and die." 

In the alternative topological scenario, it is instead the 
momentum dependence of the mass operator E(p, e) that 
gives rise to the QCP at a critical density where the 
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sum 1 + (dH(p, e)/9e°) changes sign. Let us recall that 
on the disordered side of the QCP, system properties obey 
standard FL theory. This fact, as we will see, is incom- 
patible with the collective scenario. Indeed, divergence 
of the effective mass imposes some restriction on the first 
harmonic /i entering Eq. (|6|). By way of illustration, 
we address a familiar 3D case in which simple manipu- 
lations performed on Eq. © at the critical QCP density 
Poo = kL/3tt 2 yield 

vf(Poo) _ M 1 o 

-^r = W( P ^j- 1 -3 FAp °° ) - > (9) 

where v F — pp/M. 

In principle, given the bare interaction potential V, 
the harmonic /i can be evaluated within a microscopic 
approach to FL theory (for examples, see Refs. Il6ll5ll ). 
in which this quantity is a continuous functions of the 
density p. Then the left-hand side of Eq. ([9]) must 
change sign beyond the point where the relation © is 
met, implying that the Fermi velocity vf, evaluated with 
the Landau quasiparticle momentum distribution uf(p), 
necessarily changes sign at the QCP. This behavior con- 
flicts with the collective scenario for the QCP, in which 
such a sign change is impossible. We conclude that if FL 
theory is applicable on the disordered side of the QCP, 
the collective scenario is not relevant to the QCP. The 
sign change results in the violation of the necessary sta- 
bility condition for the Landau state. Nevertheless, as 
we will see in Sec.V.B, the original Landau quasiparticle 
picture holds on both sides of the QCP, implying that in 
contrast to the collective scenario, the QCP itself and the 
line T/v(/o) of putatively associated second-order phase 
transitions are separated from each other. 

In some influential theoretical articles devoted to the 
physics of the QCP, it has been claimed that switch- 
ing on the interactions between particles fails to pro- 
duce a significant momentum dependence in the effec- 
tive interaction function /, and hence that the topologi- 
cal mechanism is irrelevant. This assertion cannot with- 
stand scrutiny. The natural measure of the strength of 
momentum-dependent forces in the medium is provided 
by the dimensionless first harmonic of the interaction 
function /(p, pi) of Landau theory. To avoid any misun- 
derstandings it is worth replacing the divergent effective 
mass M* appearing in the dimensionless Landau har- 
monic Fi = fippM* /it 2 by the free-particle mass to ob- 
tain F® = fippM/w 2 as a relevant characteristics of the 
strength of the momentum dependence of the effective in- 
teractions. In a system such as 3D liquid 3 He where the 
correlations are of moderate strength, the result F® > 2.0 
for this measure extracted from specific-heat data is al- 
ready rather large. The data on 2D liquid 3 He are yet 
more damaging to the claim of minimal momentum de- 
pendence, since the effective mass is found to diverge in 
dense films ^r— implying that F® > 3.0. In consider- 
ing the occurrence of QCP phenomena in strongly cor- 
related systems of ionic crystals, it should be borne in 
mind that the electron effective mass is greatly enhanced 



due to electron-phonon interactions that subserve po- 
laron effects^— 

The irrelevance of the collective scenario may be drawn 
with respect to the notion that a Pomeranchuk instabil- 
ity drives the QCP. To be definite, the argument may 
be framed in terms of a ferromagnetic phase transition, 
for which the associated Pomeranchuk stability condition 
reads 



I + Go(poo) = 1 + go(Poo) PooM [ iPco) = 0. (10) 



Here we set p = p^ because in the collective scenario, the 
point of divergence of the effective mass coincides with 
the point of a second-order phase transition. 

The crucial point is that at the QCP, both the rela- 
tions and (|10j) must be met simultaneously. Thus if at 
the QCP, the Pomeranchuk stability condition is violated 
then the right-hand sides of both Eq. (fTUf and Eq. (|9]) 
must vanish at the same density poo ■ However, this van- 
ishing can only happen accidentally or with the aid of 
an external magnetic field. Otherwise, we encounter a 
dilemma: either ^(/Ooo) = 3 and hence the effective 
mass diverges while the spin susceptibility remains finite; 
or else 1 + G (p oo ) = 0, and the spin susceptibility di- 
verges while the effective mass remains finite. Anyway, 
in both the cases, the end point of the line Tn(p) is sep- 
arated from the QCP. Such a separation of phases was 
first observed in 2D liquid 3 Hei2 and was also recently 
documented in precision experiments on heavy-fermion 
metals^^ 

In studies of the class of heavy-fermion metals, the in- 
formation most valuable for an understanding of their 
properties is drawn from studies of those compounds in 
which localized magnetic moments of ions forming the 
crystal lattice interact with narrow electron bands lo- 
cated near the Fermi level. The full theory of this so- 
called Kondo lattice problem is far from complete (for ex- 
ample, see the monograph of Hewson. 57 ) During the last 
decade certain results from this theory were applied in 
the analysis of the quantum critical phenomena in heavy- 
fermion metals i 15 ' 58 " — In the corresponding scenario for 
the QCP, called the Kondo breakdown model, the Kondo 
resonance is destroyed at the QCP, causing the Kondo 
temperature Tk to vanish and the FL picture to be recov- 
ered on the disordered side of the QCP. Quantitatively, 
this model, with the dynamical exponent Z = 3j£i differs 
from the standard spin-density wave (SDW) model where 
Z = 2, providing more realistic critical indexes. Notwith- 
standing this fact, the Kondo breakdown scenario for the 
QCP is burdened by the same shortcomings as the SDW 
scenario. In all, the weight of the above considerations 
compels us to "place the cart and the horse" in develop- 
ing a viable theory of the QCP, that is, to shift the focus 
from collective mechanisms to topological rearrangement 
of the Fermi surface as the origin of NFL behavior and 
the key player in the QCP regime. 
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V. EMERGENCE OF FLAT BANDS AS A 
UNIVERSAL FEATURE OF STRONGLY 
CORRELATED FERMI SYSTEMS 

In standard FL theory, applicable to weakly and mod- 
erately correlated Fermi systems, the T = ground-state 
momentum quasiparticle distribution n(p) is the Fermi 
step tif(p) — 9{pf — p)- In this case, solution of the 
Landau relation (J9j) for the quasiparticle spectrum e(p) 
is obviated, being reduced to a simple momentum inte- 
gration. The solutions for e(p) prove to be continuous 
functions, with all derivatives existing everywhere in mo- 
mentum space. On the other hand, in strongly correlated 
Fermi systems exhibiting a QCP, the basic FL formula ([T]) 
becomes fallacious, together with the asserted property 
n(p) — np(p). As we shall see, the necessary condition 
for stability of the Landau state is inevitably violated at 
some point in the QCP domain, implying that the min- 
imum of the functional E [n] is displaced from the Fermi 
step np{p) to another point in the admissible space of 
momentum distributions specified by the Pauli restric- 
tion < n(p) < 1. 



A. Case study: Rearrangement of the Landau state 
in a simple model 

In practice, numerical solution of the basic relation ([5]) 
should demonstrate whether the true T = ground-state 
momentum distribution n*(p) coincides with the Fermi 
step uf{p) or with something more interesting. How- 
ever, we need not resort to cumbersome numerical cal- 
culations to elucidate the rearrangement of the Landau 
state in strongly correlated Fermi systems. To illustrate 
this point, consider as an example a simple model where 
the energy functional E{n) has the form 



E(n) = Tr J 



P 

2M 



i{p)dv 



^TnTra j n( Pl )n(p 2 ) f ( Pl - p 2 )d Vl dv 2 (11) 



that results in the following relation 



e(p) 



P 

2M 



f i + J f(p-p 1 )n{p 1 )dv 1 (12) 

between the quasiparticle momentum distribution n(p) 
and the single-particle spectrum e(p). In what follows we 
focus on key factors that promote a rearrangement of the 
Landau state as the correlations gain strength, assuming 
the interaction function to be as follows: f(q) = X/q, 
with the zero and first harmonics coinciding with each 
other: /o = fi = X/pf. It is this form of / that was 
employed in the first article^ where the phenomenon of 
flattening of single-particle spectra was uncovered. Ac- 
cording to Ref. [0, the necessary stability condition for 
the Landau state is violated at the critical value 
f lPF M _ A/Aqcp _ „ 



F?( Po 



(13) 



Comparison of this result with Eq. (|9]) shows that the 
effective mass M* does diverge at this point. 

A straightforward estimate for the FL ground-state en- 
ergy per particle is obtained by inserting the Fermi step 
rip into Eq. (ITT1) to yield 



E(uf) 3p 



N 



WM 



1 



2F? 



(14) 



If F® is large, the interaction term is dominant. To di- 
minish this term and thereby lower E(n), the momentum 
distribution must be spread as much as possible, until 
this is rendered counterproductive due to the consequent 
increase of the kinetic term. 

To estimate the effect of a rearrangement from n(p) 
to n*(p), let us assume the new ground-state momen- 
tum distribution takes the form n*(p) = v9(pf — p), 
with v < 1. Particle- number conservation then implies 
Pf = u~ x /^pf, and the kinetic energy is increased by 
a factor (pf/pp) 2 . If the interaction function / were 
momentum- independent, no rearrangement would occur. 
In the model posed, however, f(q) falls off as 1/q, and 
therefore the replacement — ¥ n*(p) leads to the 
following result 



E(r 



3p 



N 



1QM 



El + ^EL F o 

P F 3 P/ 1 



(15) 



This energy estimate is minimized with respect to the 
parameter pf by the ratio 



EL 

Pf 



Fl 
3 



1/3 



(16) 



thus fixing v. This ratio must exceed unity, lest the 
Pauli restriction be violated; hence we surmise that the 
results obtained are meaningful provided F° > 3. Re- 
ferring to Eq. (|9]), we then infer that if this inequal- 
ity is met, the Fermi velocity vp does change its sign, 
and the Landau state is indeed destabilized. With these 
results, the energy advantage of the rearranged state, 
AE* = E(n#) — E(np), is readily estimated: 



AE* 



3p| 
10M 



Fl 
3 



1/3 



Fl 
3 



1/3- 



(17) 

Since the variational procedure is performed on a wider 
class of solutions, the rearrangement ensures a negative 
supplement AE*. 

Naturally, one cannot clarify all details of the rear- 
rangement of the Landau state by means of simple varia- 
tional estimates. It will be seen below that in correlated 
Fermi systems where the interaction function is strongly 
repulsive, the momentum distribution n* that minimizes 
the functional E(n) does in fact dive inside the defini- 
tional domain D of n(p), rather than switching between 
one boundary point of this domain and the other. In this 
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FIG. 1: Quasiparticle filling n(p) calculated for k = 0.07jjf 
and different values of coupling constant A in the model with 
interaction /(g) = \(q 2 + re 2 ) -1 . 



case, the true ground state momentum distribution n*(p) 
is found as a solution of the variational equation^ 



5E{n(p)) 
6n(p) 



-fJL = 0. 



(18) 



As a concrete illustration, Fig. [TJ presents results from 
numerical calculations performed in Ref. [l^ for the quasi- 
particle filling determined from Eq. (|18[) , employing 
f(q) = \(q 2 + k 2 )- 1 with k = 0.07p F . The resulting 
quasiparticle momentum distribution n*(p), smooth over 
the whole region < p < pf , is found by means of a regu- 
larization procedure that is effective in treating ill-posed 
problems. As applied here, this procedure involves the 
introduction of a very small additional local interaction 
between quasiparticles of the form Uo5(p — pi). Eq. (fT8|) 
then has the specific expression 



JL_ + fWp) + A |_2^ 



ra*(pi) 



dvi-fj, = 0, (19) 



) 2 + K- 



and was solved by a conventional grid method, taking 
Uo/e a F — 10~ 3 . As seen in Fig. [TJ a smooth solution of 
Eq. (|19p exists for every A value. However, this solution 
does not satisfy the Pauli principle until A attains a crit- 
ical value A cr = 1.52 Aqcp, where Aqcp corresponds to 
the first topological transition associated with the QCP. 
At A > A cr , the Pauli principle becomes immaterial. 
In the illustrative case considered, the Fermi volume is 
found to expand or swell by as much as 30%. 



B. Flattening of the spectrum of single-particle 
excitations in strongly correlated Fermi systems 

One striking feature of the phenomenon of rearrange- 
ment of the Landau state at A > A cr is a profound flat- 
tening of the spectrum of single-particle excitations e(p). 
The emergence of the flat band can be illuminated with- 
out elaborate calculations. The flattening of e(p) stems 



from the necessary condition for stability of the ground 
stated which requires that the change of the ground- 
state energy Eq produced by any admissible variation of 
the momentum distribution n(p) be nonnegative 



SEq = / e(p; n,F)5n(p)dv > 0. 



(20) 



Starting from the Landau state, the admissible variations 
of n(p) from uf{p) must obey the Pauli principle, be- 
ing negative at p < pp and positive at p > pf- In the 
standard FL theory, the signs of e(p) and p — pf always 
coincide and therefore the Landau state is stable. How- 
ever, when we address the rearranged filling of momen- 
tum states represented by the quasiparticle occu- 
pancies differ from the FL values and 1 in an interval 
[pi,Pf]] accordingly the admissible variations of n„(p) in 
this range do not have a definite sign. The only way to 
avoid violation of the inequality (|20[) in this situation is 
to make the spectrum e{p) completely flat in the domain 
Pi < p < Pf . This phenomenon, originally called fermion 
condensation^ and later swelling the Fermi surface^ or 
emergence of flat bands, was discovered 20 years ago. It 
was recently rediscovered by Lee^ in a different context, 
investigating the finite-charge-density sector of conformal 
field theory (CFT) based on the AdS/CFT gravity /gauge 
duality. 



C. Breakdown of the analyticity of solutions of the 
Landau equation for the spectrum e(p) 

From the mathematical standpoint, the FL functional 
^fl(^) is equivalent to any functional of the Thomas- 
Fermi (TF) family of functionals that depend on the 
density p(r), even though E-p^ri) and Etf(p) refer to 
different definitional domains of the arguments n and p. 
For the ensuing development, it is convenient to write 
Etf(p) in the generalized form 

E TF (p) = T {p{r)) + ( p(r)U(r)dr 



p(r 1 )p(r 2 )y(r 1 -r 2 )dridr 2) (21) 



with the normalization condition 



J p{r)dv = Z, 



(22) 



where the finite number Z is the atomic charge. In 
Eq. (|21l) . we identify r(p) as the electron kinetic energy, 
U(r) as an external field, and V(r) as the two-particle 
interaction, vanishing for r — > oo. The minimum of the 
functional (|21|) is determined from the standard varia- 
tional condition 



_ SE TF (p) 
Sp(r) 



(23) 
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The solution p(r) of this equation must be everywhere 
nonnegative; in addition, it must vanish at r — > oo to 
meet the normalization condition (f22j) . Explicitly, one 
finds 



2M 



p( r i) 
r - n| 



jdri — fi = 0, 



(24) 



where pjr(r) = (37r 2 p(r)) 1 / 3 . 

The density p(r) is known to decay rapidly as r in- 
creases beyond the TF radius ~ l/(Me 2 ), implying that 
the integral on the left side of Eq. (|24|) and the external 
field cancel each other at distances greatly exceeding this 
radius. Since p vanishes in a neutral atom, the asymp- 
totic conditions are met, and there exists a solution p(r) 
that is everywhere continuous, with all its derivatives 
defined^ 

The preceding argument breaks down for the case of 
an ionized atom. This is evident already from the singly- 
ionized case: all terms on the left in Eq. (|2~4"|) vanish 
as r — > oo but one, namely the chemical potential p, 
which has a finite valued Thus, Eq. (|24p cannot be sat- 
isfied at r — ^ oo. We conclude that this equation only 
applies within a finite domain defined by r < R, with 
R treated as the boundary of the ionized atom beyond 
which p(r) vanishes identically.— The equations deter- 
mining the minimum of the functional (|21[) now become 



P 2 F (r) 
2M 



Z-+e 2 

r .1 \r — Ti\ 



p(r) = 0, r > R. 



dri — p — 0, r < R, 

(25) 



Having obtained the corresponding density distribution, 
the function e(r) of Eq. |23|) is evaluated in closed form 
as 



e(r) = , r<R, 



e(r) = -Z^--+e 2 / ,_f v x _' - p, r > (26) 



r-ri 



The energy e(r) vanishes identically only at r < R, 
whereas at r > R it changes somewhat with the obser- 
vation point r. Since any analytic function that vanishes 
identically in some domain must be identically zero ev- 
erywhere, we see that analytic solutions of the TF vari- 
ational problem simply do not exist in the case of an 
ionized atom. 

In the preceding development, we have demonstrated 
that the Thomas-Fermi problem is ill posed, by witness- 
ing the breakdown of analyticity due to the impossibility 
of reconciling (i) the variational equation that is to de- 
termine the solution by functional minimization with (ii) 
the prescribed asymptotic behavior of the desired solu- 
tion. Similarly, the task of finding a satisfactory min- 
imum of the Landau functional £fl(?i) also belongs to 
the class of ill-posed problems. Since the first term in 



Eq. ([19]) diverges at p — > oo, there are no solutions of 
the Landau variational problem that are continuous with 
derivatives everywhere in momentum space. Therefore 
the domain of solution must be restricted, as in the TF 
problem (though in momentum space rather than coor- 
dinate space) leading to the following equation analogous 
to first of Eqs. 



|U + //(P- Pi)n*(pi)dvi - p = 0, peC. (27) 

A NFL solution n*(p) continuous in the momentum do- 
main C shows no Migdal jump at the Fermi momen- 
tum p = pp. Correspondingly, in this domain the ba- 
sic property n 2 (p) = n(p) obeyed by the standard T = 
FL quasiparticle momentum distribution no longer holds. 
Outside the domain C, the spectrum e(p) ceases to be flat, 
and accordingly the distribution n* (p) coincides with the 
Fermi step, being 1 for occupied states and 0, for empty 
ones. As stated above, the phenomenon of "swelling" of 
the Fermi surface implies that it expands from a line to 
a surface in 2D, and from a surface to a volume in 3D. 
In other words, beyond the point of fermion condensa- 
tion, the single-particle spectrum e(p) acquires a com- 
pletely flat portion e(p) = in the momentum interval 
Pi < p < pf (the so-called fermion condensate (FC)). 

With the aid of the variational condition (fT5|) . one 
can understand that the system with the FC cannot be 
treated as a gas of interacting quasiparticles, because in 
case a particle with momentum p E C is added to the 
system, the addition results in a rearrangement of the 
whole distribution function n* (p) in the C region. This is 
an essential aspect of the modification of the FL quasi- 
particle picture required to describe strongly correlated 
Fermi systems. 



D. Systems hosting a FC as a special class of Fermi 
liquids 

The flattening of the spectrum e(p) in systems with 
a FC is reflected in the structure of the single-particle 
Green function G(p, s) on the imaginary axes in the limit 
e — » 0, where^ 



G(p,e) = ~, peC. 

£ 



(28) 



This makes the essential difference in evaluation of the 
topological charge, given in the 2D case by the integral 2 -^ 

N = J {Gip^e))- 1 diGfae) — , P eC, (29) 



where the Green function is considered on the imaginary 
energy axis e — i£ and the integration is performed over a 
small contour 7 in the (p,£)-space embracing the Fermi 
line. For conventional Fermi liquids and states with a 
multi-connected Fermi surface, the topological charge Af 
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is an integer, whereas for a states having a completely flat 
portion in the spectrum e(p), integration over the contour 
embracing the Fermi line yields a half- odd-integral value 
of N£ Accordingly, from the topological point of view, 
systems hosting a FC form a special class of Fermi liquids. 

In dealing with the whole Fermi-surface region, 
Eq. (f2"5|) needs a slight generalization. To explain this 
modification, we introduce a minute temperature T. It 
is then natural to write the FC Green function in the 
form 

<-Hp,£,T) = 1 ; — , p€C, (30) 

e+«7 e— 

familiar in low-temperature physics, the occupation num- 
bers n*(p, T) being given by the Landau formula In 
writing this formula the ratio 7/e is assumed to be being 
much less than unity, which applies at least for systems 
having a small proportion of FC. 



E. Low-temperature expansion in systems with a 

FC 

It is a key feature of the phenomenon of fermion con- 
densation that a tiny elevation of temperature has prac- 
tically no effect on the T = momentum distribution, so 
that Eq. (Q} can be inverted to obtain 3 

e(p,T^0)=Tln 1 ~ n ; ( f\ P eC. (31) 
^* \P) 

This means that at extremely low temperatures, the dis- 
persion of the single-particle spectrum e(p) in the FC 
domain is proportional to T, thus allowing one to iden- 
tify (or redefine) the FC as the totality of single-particle 
states having a linear-in-T dispersion of their associated 
spectrum e(p, T). As a matter of fact, the formula (l3lT) is 
instrumental to developing low-temperature expansions 
of the basic quantities e(p,T) and n(p, T). To clarify 
how such expansions are generated, we first decompose 
the quasiparticle momentum distribution n(p,T) as fol- 
lows: 

n(p,T) = n*(p)+Tp(p,T), (32) 

where n.*(p) is the T = solution of Eq. (fT2)) . To facil- 
itate the analysis we assume that the FC occupies that 
part of momentum space lying between p — to p = pf , 
so that n*(p) vanishes identically at p > pf. 

Upon inserting the decomposition (|3"2")l into Eq. (TT2")) 
and replacing there e(p) by 7Tn[(l — n(p))/n(p)] we ar- 
rive, as before, at the system of two equations 

In 1 n ^p P ^ = J f(p,Pi)v(pi,T)dv!, p<Pf, 
-1nv(p,T) = z»(p) + f(p,p 1 )v(p 1 ,T)dv 1 , p > p/,(33) 



where = e*(p)/T, with 

p 2 f 

e*(p) = + / /(p,Pi)«*bi) rfw i ~M> P > Pf- ( 34 ) 

In deriving the second of Eqs. (|33|) . the logarithmic term 
has been simplified by the replacement 1 — v(jp) — > 1, 
noting that the quantity v{jp) is small everywhere in mo- 
mentum space. The behavior of v(p) is readily explicated 
in two regions of momentum space on either side of the 
boundary point pf. In the FC region, upon neglecting v 
on the left-hand side of the first of Eqs. (|33|) . this equation 
becomes linear, and upon introducing the low-T expan- 
sion, we have 

v(p, T) — a + a.\T + . . . , peC. (35) 

The coefficients Qo ~ 0(1), ai,... may be evaluated 
in closed form. This expansion begins to fail near the 
boundary point pf, where the ratio /n*(p)/T becomes 
less than unity and the logarithmic structure of the left- 
hand members of Eq. (|33[) exhibits its influence in the 
emergence of an additional factor ln(l/T) in v(T). Out- 
side the FC region, the solution 

z/(p,T)~e- £ * (p)/T (36) 

is easily found, provided e*(p)/T ^> 1. Near the bound- 
ary point pf where e*(p)/T ~ 1, the simple formula (|36|) 
fails, and the quantity v(T) again experiences a moderate 
enhancement, which can be evaluated only numerically. 

The impact of such enhancement on the thermody- 
namic quantities of a system hosting a FC is not signifi- 
cant, since the range of the domain where the logarithmic 
enhancement of v(T) takes place is too small. As an ex- 
ample, let us evaluate the specific heat. We recall that 
at the QCP, the Sommcrfeld ratio 

Wr^l* (37) 

diverges at T — > due to the divergence of the den- 
sity of states, which is proportional to the effective mass 
M*(T). In systems with a FC, this divergence is seen to 
disappear because the major term n* (p) of the sum (|32|) , 
being T- independent, does not contribute to 7(T). With 
this result, upon inserting Eq. (I3"5j) into Eq. (|3"T1) . the 
contribution of the FC region to j(T) is evaluated in a 
straightforward matter, retaining the first term ao — 1 of 
the expansion (I35[) . Since, as always in thermodynamic 
calculations, the condition |e|/T ~ 1 applies, the result so 
obtained turns out to be T-independent just as in stan- 
dard FL theory a characteristic feature of systems with 
a FC, first disclosed in Ref. 0. 

Let us us now estimate the contribution to the Som- 
merfeld ratio 7 from an integration region in which the 
function v is logarithmically enhanced. Since in this re- 
gion, the limits of integration are restricted by the same 
condition |e| ~ |e*| < T, we surmise that the magnitude 
of this contribution is determined by the product of the 
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factor T and the average value of the factor dp/de„(p). 
Analysis shows that the latter is indeed enhanced in the 
boundary region; one has dp/de*(p) oc T~ Q , but with 
a < 1. Collecting all factors, we find that at T — > 0, the 
contribution to 7 of the regions adjacent to the bound- 
ary point pf is suppressed as T 1- " ln(l/T). Accordingly, 
we conclude that the emergence of a FC in significant 
proportion produces a great suppression of the QCP be- 
havior of the Sommerfeld ratio. 



F. Fate of the Landau-Luttinger theorem in 
systems with a FC 

Here we demonstrate that the Landau-Luttinger the- 
orem (JTJ) breaks down in systems hosting a FC. It will 
be instructive to pinpoint the flaw that emerges in pur- 
suit of the standard Luttinger argument^ The analysis 
begins with expression of the particle number N as the 



mtcgra 



39 



N 
V 



= -2 



In 



G(p,e) 



G*( P ,e) 



dp de 

(27r) 4 i 



(38) 



In conventional Fermi liquids, this integral is propor- 
tional to the difference In (G(p,£ = 0)/G*(p,e = 0)) - 
ln(G(p, £ = — oo)/G*(p, e = —00)). In turn, the value of 
the ratio In (G(p, e)/G*(p, e)) is proportional to the argu- 
ment (f(e) of the complex function In G(p,e). One has 39 
tp{— 00) = 7T independently of the presence or absence of 
a FC in the system. Since ImG _1 (p, e) — 0, a nonzero 
value of the integral in momentum space is obtained by 
integrating over a domain where ReG(p, e = 0) > 0, 
yielding the customary result p^,/3ir 2 . However, in the 
FC region, the function ReG _1 (p, e — 0) vanishes iden- 
tically, and for states belonging to the FC, the result, 
contained in the expression 



N 
V 



pt 



Pi 
3tt 2 



d_ 

de' 



In 



G( P ,e) 
G*(p,s) 



dp de 

(2tt) 4 z 



(39) 



becomes uncertain. It is treated by dividing the region of 
energy integration into two intervals, namely from —00 
to a small negative value 5, and from <5 to 0. Accord- 
ing to Eq. (|30]). G(p,e — 5) is simply 1/e, in agree- 
ment with the property (|2"5)) . Then the integral over 
the first of the two intervals vanishes, since, at e = S, 
the argument of the ratio ln(G(p, e = S)/G*(p, e = 8)) 
coincides with that at e = —00. As a result, we are 
left with the integral over the second interval, containing 
ln(G(p, e = 0)/G*(p, £ = 0)), which can be evaluated by 
going around the pole of the Green function (|3T)]) . As a 
result, we are led to the Landau formula @, which now 
reads explicitly 



N 
V 



Pi_ 
3tt 2 



dp 



(40) 



Thus we conclude that in systems having a FC, the 
Landau-Luttinger theorem ([7]) is violated, yet the Lan- 
dau postulate ([2]) remains intact. 



VI. MERGING OF SINGLE-PARTICLE LEVELS 
IN FINITE FERMI SYSTEMS 

A phenomenon analogous to swelling of the Fermi sur- 
face in infinite Fermi systems also exists in finite sys- 
tems. The systems implicated include atomic nuclei, 
whose quantitative description within the FL framework 
was elaborated by A. B. Migdal£ In their conventional 
wisdom, textbooks teach us that under variation of input 
parameters, two neighboring single-particle levels may re- 
pel or cross one other. However, this familiar dichotomy 
overlooks a third alternative: levels can in fact merged 

This phenomenon, missing in the theory of finite Fermi 
systems, is made possible by the variation of single- 
particle energies with level occupation numbers — a prop- 
erty central to Landau theory. The primary condition 
for merging to occur is that the Landau-Migdal interac- 
tion function / is repulsive in coordinate space, which 
holds for the effective nn and pp interactions in the nu- 
clear interior— and for the electron-electron interaction in 
atoms. 



A. Schematic model of merging of single-particle 
levels 

Following Ref. 66, consider a schematic model involv- 
ing three equidistant neutron levels, separated by an en- 
ergy distance D in an open shell of a spherical nucleus. 
The levels are denoted — , 0, and +, in order of increasing 
energy. The single-particle energies e\ and wave func- 
tions v 5 a(i") = i2ni(r)$j/ m (n) are solutions of 



[p 2 /2M + E(r, p)]cpx(r) = e x <p x (r), 



(41) 



where S denotes the self-energy. In even-even spherical 
nuclei, which in their ground states have total angular 
momentum J — due to pairing correlations, the ener- 
gies e\ are independent of the magnetic quantum number 
m associated with the total single-particle angular mo- 
mentum j. We suppose that the level — is filled, the 
level + is empty, and N neutrons are added to the level 
0, changing the density p(r) by Sp(r) = N Rf lolo (r) / An . 

In what follows we shall retain only a major, spin- and 
momentum- independent part V of the self-energy E and 
a dominant, ^(r)-like portion of the Landau-Migdal inter- 
action function /. Accordingly, the FL relation between 
S and p responsible for the variation of e\(n) with n 
reduces tc£ 



SV(r) = f\p(r)]Sp(r) 



(42) 



When particles are added to the system, all energy 
levels are shifted somewhat, but the level that receives 



the particles is affected more strongly than the others. 
For the sake of simplicity, the diagonal and nondiagonal 
matrix elements of / are assigned the respective values 

u = j R 2 nl {r)f[p{r)]R? nl {r)r 2 dr/^, 

w = J R 2 nl (r)f[p(r)]R 2 nih (rydr/4n, (43) 

independently of the quantum numbers nl, ri\l\. 

Based on these assumptions and results, the dimen- 
sionless shifts £ k (N) = [e k (N] - e fc (0)] /D for k = 0, +, - 
are given by 

(o(N) = n U, (+(N) = C-(iV) = n W, (44) 

where n k — Nu/{2jk + 1) is the occupation number of 
the level k, U = u(2j Q + l)/D, and W = w{2j + 1)/D. 
It is readily verified that if fppM/w 2 ~ 1, where pp = 
\J2Mtp and tp is the Fermi energy, then the first of the 
integrals (j4"3"f has a value u ~ ep j A and therefore U ~ 1, 
since D ~ epjA 2 l z in spherical nuclei. 

According to Eqs. flU} at [U - W) > 1, the dif- 
ference d(N) = 1 + £+(N) — £,o{N) changes sign at 
noc = 1/(U — W), before filling of the level + is com- 
plete. At no > hq c , in the standard scenario provided by 
Hartree-Fock theory, all added quasiparticles must reset- 
tle into the empty single-particle level +. However, not 
all of the quasiparticles can take part in the migration 
process, since the situation would then be reversed, and 
the roles of the levels interchanged: the formerly empty 
level, lying above the formerly occupied one, would have 
the maximum positive energy shift, rendering migration 
impossible. Thus, the standard Fermi-liquid filling sce- 
nario, which prescribes that one and only one single- 
particle level lying exactly at the Fermi surface can re- 
main unfilled, while all others must be completely occu- 
pied or empty, encounters a contradiction. 

This contradiction is resolved as follows^ Migration 
occurs until the single-particle energies of the two levels 
in play coincide. As a result, both of the levels, and 
+, become partially occupied — an impossible situation 
for the standard Landau state. Solution of the problem 
reduces to finding the minimum of the relevant energy 
functional 

E a = e {0)N + e+{0)N + + ± [u(N 2 + JV*) + 2wN N + ] 

with Nh = Y] m nkrm through a variational condition 
5E 8E 

= = [i, Vm,mi, (46) 

on 0m 6n +mi 

where fj, is the chemical potential. Eqs. (|4"6")l are conve- 
niently rewritten as conditions 

e (N) = e (Q)+N o u + N + w = fi, 

e+(N) = e+ (0) + N w + N+u = fi (47) 



x 10 x 1 




x 10 x 1 



FIG. 2: Top panels: Dimensionless distance d — (e+ — £o)/D 
between levels + and as a function and as a function of 
the ratio x = 7V/(2jo + 2j+ + 2). Lower panels: Occupation 
numbers nk for levels and +. Input parameters: U = V = 
3, W = 1. For the left column, the ratio r = (2j + l)/(2j + + 
1) = 2/3; for the right, r = 3. 



for coincidence of the single-particle energies eo and e + , 
which, at N > N c = (2j + 1)/(U - W), yield N = 
±(N + N c ) and N+ = |(iV - N c ). 

Results from numerical calculations are plotted in 
Fig. [21 which consists of two columns, each made up 
of three plots. The upper panels show the dimension- 
less ratio d(x) = ]e+(#) — eo(x)] /D versus x — N/ (2 jo + 
2j + + 2) e [0, 1]. The lower panels give the occupation 
numbers n + (x) and rio(x). We observe that there are 
three different regimes: in two of them d ^ and there 
exist well-defined single-particle excitations, and in the 
third, the energies of the levels and + coincide at zero. 
Passage through the three regimes can be regarded as a 
second-order phase transition, with the occupation num- 
ber n + treated as an order parameter. 

The single-particle levels remain merged until one of 
them is completely filled. If the level fills first, as in 
the left column of Fig. [2l then under further increase 
of N, quasiparticles fill the level +, signaling that the 
distance d(N) again becomes positive. This behavior re- 
sembles the repulsion of two levels of the same symme- 
try in quantum mechanics, although here one deals with 
the single-particle levels of different symmetry. In a case 
where level + becomes fully occupied before the level 
does, as in the right column, the distance d(N) becomes 
negative, and the two levels just cross each other at this 
point. Interestingly, analogous behavior has been discov- 
ered recently^ in numerical calculations of single-particle 
spectra of multilayered systems. 
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B. Competition between shell-model effects and 
merging of single-particle levels 

As is well known that the stability of atomic nuclei 
hinges critically upon negative shell corrections to the 
Weizsacker formula for ground-state energies. These cor- 
rections are maximal in magic nuclei such as Sn 132 and 
Pb 208 , where the density of single-particle levels near the 
Fermi surface drops sharply due to the existence of an 
energy gap between the last occupied and first unoccu- 
pied levels. Such behavior is in complete agreement with 
the idea of an atomic nucleus "as a gas of interacting 
quasiparticles" £ The magnitude of shell effects decreases 
with increasing atomic numbers, primarily because the 
magnitude of the magic gap falls off rapidly in heavy and 
superheavy nuclei. With merging phenomenon coming 
into play, the situation changes. For one thing, merg- 
ing ensures a negative contribution to the ground-state 
energy, despite the fact that the density of the single- 
particle levels exhibits an opposing trend, reaching a 
maximum at the Fermi surface, rather than a minimum, 
as in the conventional shell model. 

In effect, the existence of gaps in filling prevents the 
merging phenomenon from working in full force. How- 
ever, in heavy and superheavy atomic nuclei where this 
gap becomes smaller and smaller, the FC correction to 
the Weizsacker ground state-energy grows in importance, 
creating the opportunity for level merging to be an im- 
portant playmaker in the formation of new stability is- 
lands. 



C. Merging of single-particle levels in an external 
magnetic field 

In this subsection, we employ the simplest model of 
fermion condensation, proposed in Ref. [H, to investigate 
the impact of an external magnetic field on the FC spec- 
trum and demonstrate that the presence of this field does 
not prevent the Fermi surface from swelling. The anal- 
ysis is restricted to the simplest form of the FL energy 
functional, namely 



If 1 
E= / r(r)dr + - 

2M K ' 2 



/(r-s)p 2 (r,s)drds, (48) 



in which / stands for the Landau interaction function. 
The density matrix p and kinetic-energy density t are 
respectively given by 

p(r,s) =Y,n&*{v)* i {s), r(r) = £ r^V^V^r), 

(49) 

in terms of the true quasiparticle occupation numbers rij 
of the single-particle states. Imposition of the external 
magnetic field H leads to a set of coupled equations 



(p - eA(r)) s 
2M 



^•(r)+ J /(r-s)p(r, B^sjds - E^j(v), 

(50) 



for the eigenvalues Ej and eigenfunctions tyj , where A is 
the corresponding vector potential. Specifically, we ad- 
dress 2D homogeneous electron systems, employing the 
Landau gauge in which A x — —Hy and A y = 0. Gener- 
ally, each of the equations (|50[) contains a sum of an in- 
finite number of terms involving different wave functions 
^i. This complication is normally surmounted by work- 
ing within a quasiclassical approximation, 39 such that 



p(r,s, A) 



A)(r - s), 



(51) 



where po is the well-known quasiparticle density matrix 
for homogeneous matter. However, the quasiclassical ap- 
proximation fails in systems having long-range interac- 
tion functions of the kind responsible for flattening of 
the single-particle spectrum e(p) beyond the QCP. 

On the other hand, it is just the long-range character 
of the relevant effective interactions that allows one to 
proceed while avoiding the approximation (j5"Tj) . Indeed, 
suppose the correlation radius £ specifying the behavior 
of / at large distances is in excess of the Larmor radius 
th = 1/VeH. In the sum over i implicit in Eq. (|50p . 
the diagonal term is dominant, and the initial structure 
^>(x, y) = e lp * x x{y) °f the corresponding wave functions 
is recovered. Neglecting corrections of order r 2 H /£ 2 com- 
ing from nondiagonal contributions, we are then left with 
the Landau-like equation 



1 

2M 



(p 2 y + e 2 H 2 {y- 



yo) 2 ) Xj(y) = (Ej - fnj)xj{v), 



where / is the effective coupling constant. As usual, 
yo specifies the location of the Larmor circle, while the 
occupation number nj — [1 + exp(e.,-/T)] _1 is itself de- 
pendent upon the deviation ej = Ej — p of the single- 
particle level from the chemical potential p. Now, ac- 
cording to Eq. (I52|) . the initial equidistant spectrum 
E° = Ej(H = 0) = cj°(j + 1/2), with w°. = eH/M, 
gives way to a modified level scheme determined by inte- 
ger solutions of nonlinear algebraic equation 



e(x) = uj®x — p + fn(x). 



(53) 



Recalling that n(e = 0) = 1/2, Eq. (|53[) may be recast in 
the form 



(54) 



where y = x — (p— f /2) juP c gives the displacement of the 
single-particle levels from the Fermi surface. Eq. (|5~4"|) is 
reminiscent of the corresponding equation for the NFL 
occupation numbers in the Nozieres model^ developed 
for homogeneous matter with an interaction function 
/(r — s) = const, of infinite radius. The QCP of the 
model 3 is located at / = 0, with the flattening phe- 
nomenon in effect at / > 0. 

Eq. (1541 can easily be solved graphically at T = (see 
Fig. This is done by plotting the kink n(e), together 
with the set of straight lines representing the right side 
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FIG. 3: Graphical illustration of the solution of Eq. (|54l) 

of Eq. (|54p. against the input parameter y. At y > r, the 
right side of Eq. (|5~4"|) has the form e — a, with a > 0. In 
this case, crossing points lie in the right half-plane e > 
where n — 0, implying that the initial Landau spectrum 
E® = cj c (j + 1/2) is unaffected by the interactions. On the 
other hand, at y < — r, the crossing points transfer to the 
left half-plane e < where n = 1, so that the equidistant 
Landau spectrum is simply shifted by the constant /. 

What happens in the interval — r < y < r? In this 
case, as seen from Fig. [31 the crossing points land on the 
vertical segment of the kink, yielding the solution e = 0; 
Eq. ([54]) reduces to relation 

n*(y) = 7j (l - ~) ) -r<y<r, (55) 

for finding the NFL occupation numbers n*(y,H) at 
H =/= 0. Comparison with the corresponding result 
n*(v) = (m — P 2 /(2A/))// for the Nozieres model 3 shows 
that the swelling of the Fermi surface inherent in this 
model persists, the range of the FC domain being un- 
affected by imposition of the magnetic field. At finite 
T, analysis of Eq. (|54[) confirms that in the FC domain, 
the derivative de(x)/dx becomes proportional to T, again 
in accordance with the results obtained in homogeneous 
matter at H = Oi 3 - 



VII. LIFSHITZ PHASE DIAGRAM: 
EVOLUTION FROM THE LANDAU STATE TO 
THE FC STATE 

Low-temperature phase transitions occurring in Fermi 
systems have already been studied for a hundred years. 
As a rule, these are second-order phase transitions, epit- 
omized in the phenomenon of superconductivity discov- 
ered by Kamerlingh-Onnes in 1911. By contrast, experi- 
mental studies of topological phase transitions (TPTs) 
have gained prominence only during the last decade. 
This disparity stems from the fact that low-temperature, 
second-order phase transitions associated with violation 
of one or another of the Pomeranchuk stability condi- 
tions may arise in weakly or moderately correlated sys- 
tems, as in the example of superconductivity. In the case 



of TPTs, it is instead the necessary stability condition 
(|2"U| for the Landau state that suffers a breakdown, an 
eventuality that can be realized only in strongly corre- 
lated Fermi liquids, notably in dense films of liquid 3 He 
and the dilute 2D electron liquid in MOSFETs. Nothing 
was known about such systems prior to recent experi- 
mental and technological advances. To our knowledge, 
the possibility of a many-fermion ground state with a 
multi-connected Fermi surface was first entertained by 
Frolich in 1950,-— Theoretical consideration of topologi- 
cal aspects of the electron liquid in solids dates back to 
an article by I. M. Lifshitz published ten years later^ 

The Lifshitz analysis requires specific information on 
the structure of the single-particle spectrum e(p), which 
is often not available because this spectrum may be al- 
tered dramatically upon switching on the effective inter- 
actions between quasiparticles. 1 For many years, insights 
into the effects of these interactions were pursued within 
the Hartree-Fock frameworki22r— It was not until 1990 
that such limited descriptions began to be replaced by 
more powerful FL approaches^"— ^ ~ 19 i 24 ~ 27 i 74 

In dealing with homogeneous systems, any TPT is 
associated with a change in the number of roots of 
equation^ 

e(p;n,p) = 0, (56) 

where e(p; n) is the spectrum of single-particle excita- 
tions measured from the Fermi surface and evaluated self- 
consistently with the T = ground-state quasiparticle 
momentum distribution n{p). In conventional homoge- 
neous Fermi liquids such as 3D liquid 3 He, Eq. (f5T>)) has 
a single root, the Fermi momentum pp', hence the cor- 
responding ground-state quasiparticle momentum distri- 
bution is simply Uf{p) = 0{pf —p)- 

The QCP corresponds to the particular case of a topo- 
logical transition for which the bifurcation point pb co- 
incides with the Fermi momentum pp. By its definition, 
at the QCP the effective mass diverges and therefore the 
Fermi velocity vf(Poo) — {de(p, Poo)/dp) vanishes. In 
this situation, the spectrum is constrained to have an 
inflection pointj22 i.e. two additional roots of Eq. (|5"6"]) 
emerge simultaneously at the Fermi surface, implying 
that e(p) tx (p — pf) 3 at the QCP itself. It will be seen 
that such a situation is not typical. As a matter of fact, as 
the system moves into the quantum critical regime from 
the FL domain, Eq. (|56[) must already acquire two new 
roots pi,P2 7^ Pf at some critical value Aqcp of the in- 
put parameter A. Accordingly, the Fermi surface becomes 
multi-connected, but with the ground-state momentum 
distribution n(p) — 6(jp\ — p) — 6{p2 — p) + 0(pf — p), 
still satisfying the usual FL relation n 2 (p) = n(p). On 
the other hand, those solutions of Eq. © emerging at 
some critical coupling Afc and featuring a completely 
flat portion of the spectrum e(p) correspond to solutions 
of Eq. (|56|) having an infinite number of roots, or cquiv- 
alently, solutions with an infinite number of sheets of the 
Fermi surface. We reiterate that in the associated FC 
phase, the fundamental relation n 2 (p) = n(p) of stan- 
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dard FL theory no longer holds; partial occupation of 
single-particle states must occur. 

Numerous calculations have demonstrated that as a 
rule, the ratio Afc/Aqcp ^ es between 1 and 2. When 
the coupling constant A increases beyond Aqcp into the 
regime of instability of the Landau state, a striking se- 
quence of rearrangements of the single-particle degrees of 
freedom is witnessed^ as pockets of the Fermi surface 
breed and proliferate rapidly, in resemblance to a frac- 
tal cascade. Our first task in the present section is to 
examine the details of the evolution of the topology of 
the ground state in the Lifshitz phase diagram over the 
coupling interval Aqcp < A < Afc ■ Our second goal is to 
expose the intricate evolution, with increasing tempera- 
ture, of the single-particle spectrum e(p) in phases with 
having a multi-connected Fermi surface. 



Generic features of topological transitions in 
Fermi liquids 



In principle, the bifurcation pt of Eq. (|56|) giving rise 
to the first TPT can emerge at any point of momentum 
space. Leaving aside the case addressed in the previous 
section, in which the transition leads to a swollen Fermi 
surface, we are left with the following options. The first 
and common option is pt = 0. The second and partic- 
ular one is the QCP for which pb = pp. In the third 
case < pb < pf, the bifurcation point pb is a local 
maximum of the spectrum e(p), since one has e(p) < 
throughout the Fermi volume. This implies that between 
Pb and pf there also exists a local minimum of e(p), be- 
cause by definition the spectrum e(p) changes sign at 
p = pf- Thus, if the bifurcation point is located be- 
tween and pf, the equation de(p)/dp — has at least 
two roots in this momentum interval. Simple analysis 
demonstrates that existence of at least two roots of this 
equation is also inherent in the fourth and last case, in 
which Eq. (|56l) bifurcates at a momentum greater than 
Pf and the bifurcation point is a local minimum of the 
curve e(p). 

A prerequisite for the occurrence of the bifurcations 
in Eq. (|56p is a marked momentum dependence of the 
interaction function /(p, Pi). Such a momentum depen- 
dence is often associated with forces of long range in in 
coordinate space that are induced in advance of an im- 
pending second-order phase transition. As a pertinent 
illustration, let us consider the situation in the vicinity 
of pion condensation, a phase transition in dense neu- 
tron matter predicted by A. B. Migda h 75 ' 76 We restrict 
ourselves to the case of moderate pion fluctuations where 
the Ornstein-Zernike approximation, allowing convenient 
treatment of the fluctuation contribution to the interac- 
tion function /, is still applicable. In this case, the ex- 
change contribution f n7r to the interaction function / has 
the formS 
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*Hp) + (q 2 /q 2 - i) 



(57) 



Q. 



> 0.2 




FIG. 4: Upper panel: Position of the bifurcation point pb in 
units of the Fermi momentum, versus the critical wave number 
q c (also in units of Pf)- Lower panel: critical parameters Kb 
and functions of q c /pF- 



with coupling constant g and the critical momentum 
taken as q c ~ (0.7 -j- 1.0) pf- The stiffness coefficient 
K 2 (p) vanishes at the critical density p Cl specifying the 
point where pion condensation sets in. For simplicity, we 
choose g — l/(2m^), corresponding to bare ttNN ver- 
tices, and study the behavior of quantities versus k(p). 

In the Landau state, corrections to the neutron spec- 
trum e(p) associated with the exchange of pion fluctua- 
tions are evaluated in closed form by means of Eq. ([6]), 
based on the quasiparticle distribution np{p). Substitut- 
ing the expression (|57p and integrating, we obtain 



de(p) 
dp 



P_ 

M 



ml 



16ir 2 p 2 



1 



-Pf- 



where 



and 



L{p) = In 



[(p+PF) 2 -q 2 c f 



? 4 
« 1c 



A(p) 



[(p-p F ) 



„212 



? A 



(58) 



(59) 



A(p) = arctan 



(PF+p) 2 -q 2 

K<g 



arctan 



{PF-pf - q 2 c 



mi 



(60) 

Further integration can be performed analytically^ but 
the results are too cumbersome to present here. 

Results of numerical calculations displayed in the up- 
per panel of Fig. |4] demonstrate that a new root of 
Eq. (|5oT) . lying at the origin, appears already at k = 
Kb ~ 0.4, signaling that the Landau state has become 
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FIG. 5: Neutron quasiparticle spectra e(p) (in units of £%) 
evaluated for q c = 0.9 pf (panel (a)), q c — 1.05 pf (panel 
(b)), q c = 1.14pF (panel (c)), and q c = 1.2pp (panel (d)). 
Corresponding values of the parameter k are indicated near 
the curves. 



unstable well before the system attains the point of pion 
condensation. At customary values of the critical mo- 
mentum q c , one has pb = 0. However, as q c increases to 
greater values, pb rapidly moves toward the Fermi mo- 
mentum and crosses the Fermi surface at q c ~ l.lApp. 
The dependence of the critical parameter Kb on the wave 
number q c , depicted in the lower panel, shows that the 
largest values of Kb are achieved just in the preferred 
range q c /pF ~ 0.7— 1.0. The value k^ of k at which the 
border of the instability region reaches the Fermi momen- 
tum pf is also plotted in the lower panel of Fig. @] The 
resulting curve lies below the curve of Kb(q c ) everywhere 
except for the point of contact at q c ~ 1.14pp. 

The corresponding evolution of the neutron spectra is 
illustrated by panels (b)-(d) of Fig. [5] where the spectra 
e(p) calculated for q c — 1.05pf, L14p^, and 1.2 pp are 
drawn. We see that at q c < pp the bifurcation point pb 
does emerge quite close to the origin. This topological 
transition plays important role in the acceleration of neu- 
trino cooling of neutron starsj^l since the presence of a 
new pocket of the Fermi surface in the low-momentum 
region lifts the ban on neutron /3-decay and activates the 
highly efficient direct-Urea process. 

As seen from panels (d) and (e) of Fig. [5j beyond the 
Fermi liquid QCP the curve e(p) crosses the Fermi level 
three times, the T — Landau momentum distribution 
np(p) — 0(pf — p) is rearranged, and the Fermi surface 
becomes multi-connected. In this instance, where there 
is a single bubble in the filling pattern, the new T = 
occupation numbers n(p, T — 0) are given by an alternat- 
ing sequence of two numbers, 1 and 0. The occupancy 



is n(p) = 1 at p < pi, while at p\ < p < p 2 , there is a 
gap in filling where n(p) = 0, and at p 2 < P < P3, the 
occupancy is once again n{p) — 1. This is a typical topo- 
logical phase transition, at which no symmetry inherent 
in the ground state is violated. 



B. Prerequisites for the proximity of the first TPT 
and the QCP 

Interest to the mechanisms underlying violation of FL 
theory has been sparked by the discovery of NFL be- 
havior of electron systems present in numerous heavy- 
fermion metals. It is a specific and significant feature of 
these NFL phenomena that the bifurcation point pb lies 
close to the Fermi momentum pp. To clarify the con- 
ditions necessary for the proximity of pb and pp, it is 
natural to begin with the Taylor expansion of the group 
velocity, 



f x , P-Pf 1 (p~p F ) 
v{p) = v F + vi h -v 2 t • 

Poo 2 pp 



(61) 



applicable in the immediate vicinity of the Fermi surface, 
which yields the spectrum 

e(p) = v F {p-p F ) + 77— vi{p~Pf) 2 + ^-2-v 2 (p- pp) 3 - 
2p F 6p F 

(62) 

Upon inserting this expression into Eq. (I56|) . we observe 
(once again) that already on the disordered side of the 
QCP where vf(p) is still positive, this equation picks up 
two additional real roots: 



Pi, 2 - Pf 



-pp 



3v\ 
2v~2 



1 ± 



8vf(p)v 2 



3vl 



(63) 



Accordingly, the first TPT occurs at a critical density p t 
where the Fermi velocity is given by vp(p t ) = 3vf/8v 2 . 
This analysis is self-consistent as long as the ratio vi/v 2 
is small; otherwise at least one of the roots lies too far 
from the Fermi momentum for the Taylor expansion (|6ip 
to be valid. In conventional Fermi liquids, where /(p, Pi) 
is a smooth function in momentum space and hence the 
coefficients v\ and v 2 are of the same order, the stated 
condition is difficult to satisfy. However, the situation 
becomes more favorable in the vicinity of a second-order 
phase transition where the fluctuations are of long wave- 
length. Assuming these fluctuations to be rather weak, 
their contribution J 7 to the interaction function / can be 
accounted for within well-known Ornstein-Zernike (OZ) 
approximation to yield 



•^aVy^PiPi) g 2 < T ap<r~ (S x(\p - 



Pi I 



(64) 



Here g is the spin-fluctuation vertex and x(l) = 4ir/(q 2 + 
£~ 2 ) is the spin susceptibility, £; being the correlation ra- 
dius, which diverges at the second-order phase transition 
point. If the spin-fluctuation term is well-pronounced, 
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then £ substantially exceeds the distance between par- 
ticles. The applicability of the OZ approximation re- 
quires that the inequality (g 2 /ttv f )\x\(pfO < 1 should 
be meti^i 

The spin- fluctuation contribution v si (p) to the group 
velocity v(p) is evaluated straightforwardly, employing 



the identity a a p(T^s — ^S a s5. 



y/3 — \<y a &<y 1 p to obtain 

e(p-p F ? 



2tt 



(65) 



Then instead of Eqs. (|61|) and (|62|) one obtains (with 
v F = pf/M and V2 — X s v F ) 



p F P~PF 
V ^ = W +V1 ^F~ 



3 o (p-Pf) 2 



Pi 



/ \ PF , (P~ Pf) 2 . , o (P - Pf) C ' 

where the effective mass M* is given by 



M 



A/ 



M, 



FL 



3 ^ln(2 PF C(p)), 



(66) 



(67) 



Mp L being the FL effective mass evaluated without in- 
clusion of the fluctuation contribution (which drives the 
system toward the QCP), while 



\ - P 2 r? 9 
*s — c, Pp — n 

TTVp 



HpfO 



> i. 



(68) 



Since the condition V2/V1 cx 1/A S <C 1 holds in the 
present case, the proximity of the first TPT bifurcation 
point to the Fermi momentum is established. 

In addition, the density p^ corresponding to the topo- 
logical QCP, which signals instability of the Landau 
state, is determined by condition 



M 



9 



3M* L ( Po 



r ln(2p F £ (poo)). 



(69) 



In strongly correlated Fermi systems, the ratio M/M FL 
is suppressed, implying that the stability condition (|69[) 
is violated well before the correlation radius £(p) grows to 
values sufficient to put the OZ approximation in jeopardy. 

It is worth estimating the limits of the applicability of 
the Taylor expansion (|66l) . As seen from Eq. (|65l) , the ex- 
pansion breaks down in case the difference p — pp reaches 
values of order of 1 /£. At larger distances from the Fermi 
surface, the logarithmic term in Eq. (|65|) asserts itself, in- 
troducing logarithmic corrections to the specific heat and 
other thermodynamic quantities (as shown explicitly be- 
low). 



C. Cascade of topological transitions in the QCP 
region 

The emergence of new small pockets of the Fermi 
surface is an integral feature of the QCP phenomenon, 
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FIG. 6: The phase diagram of neutron matter in (q c , k) vari- 
ables. The Landau phase (FL) occupies the white region of 
the plane. The index i indicates the number of sheets (or 
pockets) of the Fermi surface. 



whether one deals with 2D liquid 3 He, high-T c supercon- 
ductors, or strange heavy-fermion metals. A striking pe- 
culiarity of the segment of the Lifshitz phase diagram be- 
yond the first TPT point, revealed in Refs. [l7J, is the fast 
breeding of these pockets as A increases beyond Aqcp- 
This process is illustrated in Fig. |5l which shows theoret- 
ical results^ for the phase diagram of neutron matter in 
the (q c , k) plane, from calculations performed for the in- 
teraction function (f5T|) in the vicinity of the threshold of 
pion condensation. The Landau state with n(p) = uf{p) 
occupies the white region of the diagram (labeled FL in 
the figure), while the phases exhibiting new pockets of 
the Fermi surface populate the shaded portion of the 
plane, which is se par ated from the FL domain by the 
curve Kf,(q c ). Ref. |2fJ reports interesting results from a 
simple model incorporating feedback of fermion conden- 
sation on critical fluctuations. In particular, the range 
of territory on the phase diagram in which the FC wins 
the contest with the multi-pocket phase was estimated 
as 0.95 < q c /pF < 1-3 at k — 0.1. 

As seen in Figs. [7] and [5J the breeding phenomenon also 
expresses itself vigorously in a different modei^ based on 
the interaction function f(q) — g/(q 2 + k 2 ). 

The occurrence such cascades is associated with the 
condition n(p) < 1 enforced by the Pauli principle. As 
we have seen earlier, solutions of the equilibrium equa- 
tion (|18p always exist, but these solutions violate the re- 
striction n(p) < 1 until the coupling constant A attains 
a critical value Ape- It is this restriction that triggers 
a cascade of TPTs in the segment of the Lifshitz phase 
diagram defined by Aqcp < A < Afc • Indeed, in the TF 
problem where no such restriction is present, no topolog- 
ical transitions exist. Details of the transitions between 
phases with a multi-connected Fermi surface and phases 
hosting a FC remain unclear at this point. 
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FIG. 7: Cascade of topological transitions, evaluated in the 
model with f(q) = X(q 2 + k 2 ) -1 with k = 0.07pf for four 
values of A/Aqcp = 1.08, 1.11, 1.16, 1.20. 
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FIG. 8: Single-particle spectrum e(p) in units of e% — p%/2M 
together with the self-consistent quasiparticle momentum dis- 
tribution n(p), as evaluated in Ref. within a model in 
which the Landau interaction function has the form f(q) — 
X(q 2 + k; 2 ) -1 , with k — 0.07pF. Calculated for an effective 
coupling constant A, exceeding the critical QCP value Aqcp 
by 30%, this spectrum possesses nine nodes, corresponding to 
four hole pockets in momentum space. 



In this section, we discuss results of numerical calcu- 
lations of single-particle spectra e(p, T) and quasipar- 
ticle momentum distributions n(p, T) beyond the first 
TPT. 2 - As we have seen, (i) these transitions are associ- 
ated with the violation of the necessary stability con- 
dition ([20")) for the Landau state at T = 0, in which 
n(p) = 9{pf — p), and (ii) there exist different alterna- 
tives for the rearrangement of the Landau state beyond 
the transition point. The first possibility is the emer- 
gence of a small pocket of the Fermi surface; hence, the 
new quasiparticle momentum distribution n(p, T — 0) = 
0{pi—p) — 6(p2—p) + 0(pF—p) features a combination of 




FIG. 9: Single-particle spectrum e(p) (top panels) and deriva- 
tive de(p)/dp in units of v F — pf/M (bottom panels), eval- 
uated in the model (|70[) with parameters Q2 = —0.16 and 
P2 = 0.14, chosen to adequately describe results of micro- 
scopic calculations^ of 2D-electron-gas single-particle spec- 
trum at T — 0, are shown as functions of p/pf at four (line- 
coded) temperatures expressed in units of 10 -2 £%• at r s — 6.8 
(left column), r 3 — 6.9 (middle column), and r s = 7.0 (right 
column) . 



three kinks. The second possibility amounts to a swelling 
of the Fermi surface, so that the ne momentum distribu- 
tion n(p, T = 0) is continuous over a momentum interval 
[pi;P2] embracing the boundary momentum pp- In this 
momentum interval, the spectrum e(p, T) is completely 
flat at T — 0, and, according to Eq. (pll"j) . it is inclined at 
finite T with a slope proportional to T. 

On the other hand, in the state having a multi- 
connected Fermi surface, the spectrum e(p, T = 0) varies 
smoothly in the space beyond the sheets of the Fermi 
surface but oscillates rapidly in the space between them. 
The magnitude of the departure of \e(p, T = 0)| from in 
this domain, denoted below by T m , emerges as a new en- 
ergy scale of the problem, at which "melting" of the struc- 
ture associated with the well-defined multi-connected 
Fermi surface can occur. Indeed, as seen from the basic 
Landau formula ((4]), the distribution n(p, T = 0) remains 
almost unchanged as long as T < T m . However, as the 
temperature T increases, kinks in n(p) become smeared, 
and at T > T m , in the momentum domain between the 
sheets, the function n(p, T) becomes continuous and al- 
most T- independent. In this case, according to Eq. (PT|) . 
the dispersion of the single-particle spectrum e(j>) be- 
comes proportional to T. Thus at T > T m> the tempera- 
ture evolution of e(p, T) and n(p, T) becomes universal. 

To provide tangible support for this important conclu- 
sion, Figs. l9l [T"3l present the results of model calculations 
of the spectrum e (p, T ) and momentum distributions n(p) 
performed in Ref. [26| on the basis of Eq. ©. Two differ- 
ent interaction functions / are employed in these calcu- 



17 



c 0.5 

0.0 
0.01 

S o.oo 

-0.01 



0.9 



P/P F 
1.0 



1.1 



0.9 




r =7.1 



m , 10 4 




0.9 / 1.0 \ 1.1 

Pi P/P F P, 



p/p F 

1.0 



1.1 




0.9 / 1.0 \ 1.1 

Pi P ; P F P, 



1.0 ■ 



0.5 



0.0 
0.05 



s 0.00 

CO 

-0.05 

4 

2 

t Qi 

Cl 

-4 ■ 



P/P F 

0.9 1.0 



1.1 







r =7.1 : 

s 




— - 





















0.9 / 1.0\ 1.1 
Pi P'Pp P, 




0.9 / 1.0 \ 1.1 

Pi P'P F P, 



FIG. 10: Occupation numbers n(p) (top panels), single- 
particle spectrum e(p) in units of e% (middle panels), and 
ratio e(p)/T (bottom panels) for the model (|70[1 at r s = 7.1 
(left column) and r s = 7.2 (right column), exceeding the QCP 
value Too — 7.0, evaluated with the same set of the parame- 
ters as in Fig. [9] All three quantities are shown as functions 
of p/pf at different temperatures in units of e%, lower than 
the transition temperature T m = 10~ 
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The first is chosen because this function adequately 
reproduces, within the extended Landau-Migdal treat- 
ment, the results of microscopic calculations^ of the 
single-particle spectrum e(p) of the 2D electron gas at 
T = 0, and because the corresponding initial TPT essen- 
tially coincides with the QCP located at r s — r m = 7.0. 
The second choice of / is relevant to a 3D system with 
spin fluctuations and has already been considered above 
in Sec. VII.B. 

Fig. [9] displays results for the single-particle spectrum 
e(p) and the group velocity de(p)/dp, obtained on the FL 
side of the QCP for the model ([701) . It is seen that as 
the QCP is approached, the group velocity becomes a 
parabolic function of momentum p. With increasing r s 
its minimum gradually approaches the horizontal axis, 
touching down for r s — exactly at the point p — pp, 
so that e(p,T = 0) ~ (p — pf) 3 has an inflection point 



FIG. 11: Same as in Fig. HOI but at temperatures, higher than 



at the Fermi surface. At r s > r^, the group velocity vf 
evaluated with n(p) = uf{p) changes its sign, triggering 
the rearrangement of the Landau state, in accordance 
with the topological scenario for the QCP. 

We turn now to Figs. [10] and [Til which show numer- 
ical results for the spectrum e{p) and quasiparticle mo- 
mentum distribution n(p), calculated for the interaction 
function ([70]) both below and above the temperature T m . 
It should be noted that the value of T m is somewhat un- 
certain, since we are dealing here with a crossover rather 
than a well-defined phase transition. Nevertheless, com- 
parison of these two figures reveals a striking alteration of 
the structure of both e(p) and n{p) upon passage across 
T m from the lower to the higher temperature. Indeed, as 
seen from the top panel of Fig. [I0l a well-defined multi- 
connected Fermi surface, distinguished by a pronounced 
gap in filling, exists only at extremely low temperatures 
T < 10- 4 e° F , while at T ~ T m ~ 10~ 3 e^, the gap in the 
occupation numbers closes, although at T < T m , some 
T-dependence of n(p) in the vicinity of the Fermi sur- 
face is still present. On the other hand, as seen from 
the top panel of Fig. [TTJ a T-independent behavior of 
the momentum distribution n(p), existing at T > T m in 
the FC domain, holds over a wide temperature interval. 
Comparison of the third panels of Figs. [TU] and [TTJ demon- 
strates that huge variations of the ratio e(p)/T, existing 
at T < T m , completely disappear in the FC domain at 
T>T m . 

The same conclusions follow from Figs. [12] and [T3l 
which display the temperature behavior of the spectrum 
e(p) and momentum distribution n(p), evaluated for the 
interaction function (|71[) below and above T m . 
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FIG. 12: Occupation numbers n(p) (top panel) , single-particle 
spectrum e(p) in units of 10~ 3 e% (middle panel), and ratio 
e(p)/T (bottom panel), plotted versus p/pf at three line-type- 
coded temperatures in units of £%, taken below the transition 
temperature T m — 3 x 10 -3 e%. The model (|7ip is assumed 
with parameters = 0.45 and fiz = 0.07. 



FIG. 13: Same as in Fig. [12] but at temperature higher than 



VIII. UNCONVENTIONAL 
THERMODYNAMICS OF SYSTEMS WITH A FC 



It is instructive to compare the results shown in 
Figs. IT0Hl~3l with corresponding results derived for the 
singular interaction function 



/(*) = .9, 



ir 2 e~P* k /PF 
M k 



(72) 



studied in Ref. 16 as a model of FC analytically soluble 
at T = 0. (All the input parameters in the interactions 
defined in Eqs. (|70|) . (|7T|) . and (|72|) are dimensionless.) 

It is characteristic of systems with a singular interac- 
tion function (singular as k — > 0, hence of long range in 
coordinate space) that the linear-in-T dispersion of the 
single-particle spectrum symptomatic of the FC domain 
exists at any temperature T. This behavior is evident 
in the numerical results for the interaction model ([72]) 
plotted Fig. QH 

From these and other numerical and analytical studies, 
it may be concluded that specific features of the flatten- 
ing of single-particle spectra beyond the first TPT point 
are universal, differing only with respect to the scale tem- 
perature T m . 



In this section we will focus on certain unusual features 
of the thermodynamic properties of strongly correlated 
systems hosting a FC. In making concrete estimations, it 
will be assumed that the interaction function has a small, 
specific fluctuation contribution, the origin and form of 
which has been analyzed in Sec. VII. B. 



A. Logarithmic corrections to j(T) 

The standard FL behavior of j(T — > 0) = const, 
does in fact re-emerge beyond the point of a TPT, be- 
cause the density of states N(0) again becomes finite. 
However, recovery of FL behavior occurs only at ex- 
tremely low temperatures. We will see that, as the tem- 
perature increases, 7 acquires a logarithmic correction 
Sifl{T) ex. ln(l/r). Such a correction is in fact observed 
in a number of heavy-fermion metals at temperatures 
exceeding a rather low temperature Tq. Its presence is 
usually explained within the framework of the collective 
scenario for the QCP. As already discussed, in this sce- 
nario quasiparticle weight z vanishes at a second-order 
transition point. Consequently, the mass operator S(p, e) 
contains a marginal term ~ elne, giving rise to the ap- 
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FIG. 14: Single-particle spectrum e(p) in units of e% (top 
panels), occupation number n(p) (middle panels), and ratio 
e(p)/T (bottom panels) for the model ()72[) with the parame- 
ters /3* = 10 and g* = 70.0. All three quantities are shown as 
functions of p/pf at different line-type-coded temperatures in 
units of s F . 



pearance of a logarithmic correction ^7fl(T) oc ln(l/T) 
to the FL Sommerfeld ratio 7fl(T) = const. 

However, as we will be readily appreciated, breakdown 
of the quasiparticle picture is not a necessary condition 
for the emergence of logarithmic corrections to the Som- 
merfeld ratio 7(T). Such corrections can emerge in the 
QCP region due to the presence of small spin-fluctuation 
contributions to the FL interaction function. To see this, 
let us evaluate j(T) within the spin- fluctuation scenario 
introduced in Sec. VII. B. One has 



7 (T) ~ T- 



(73) 



where n(p,T) = (1 + exp (e(p,T)/T)) is the usual 
Landau quasiparticle occupation number. Referring to 
Eq. (|65|) we know already that at <C (p — Pf) <C pf, 
the group velocity v(p) takes on a logarithmic correction 
~ In ((p — Pf)/pf) ~ In (M* c/pf), implying that 



1 
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(74) 



Since the overwhelming contributions to 7(T) and other 
thermodynamic quantities come from the energy region 



where |e| ~ T, we introduce a new variable y = e/T to 
obtain the standard FL formula 



-AT.,,) -~ / y 2 n {y){l-n{y))^-dy 



(75) 



Upon inserting Eq. (|74|) into this equation, we infer that 
the logarithmic correction to 7fl = const, becomes well- 
pronounced already at a rather low temperature To = 
t f/(pfQ, with Tp = p 2 F /2M*, according to 



S7 FL (T) 3g 2 M* 



7FL 



np F 



In 



Tf 
T 



T n < T <T F 



(76) 



The T-dependence of this correction coincides with that 
obtained within the conventional scenario for the QCP. 
But importantly: in contrast to this scenario, the regime 
involved is far from the critical one where the quasiparti- 
cle picture breaks down. Similar logarithmic corrections 
emerge in the same temperature interval when evaluating 
the thermal expansion coefficient, magnetic susceptibil- 
ity, and other thermodynamic quantities. Their presence 
is irrelevant to violation of the quasiparticle picture. 

The value of the parameter (pf£,)i an d hence that of 
the governing parameter A, can be extracted from experi- 
mental data obtained on the upper and lower boundaries 
of the interval where the logarithmic behavior ([76]) is in 
effect. The analysis of these data yields 10 < pf£, < 10 2 , 
implying that 



10 < A < 10 4 



(77) 



Comparing the Sommerfeld ratios on opposite edges of 
the interval of logarithmic behavior, one can also estimate 
the value of another parameter, namely g 2 /(Trv F ). Since 
7 drops by factor 2 or so, we conclude that g /{itVp) < 

VHpfO <o.i. 

These formulas, in first turn, Eq. (|75p . are valid pro- 
vided the inequality To < T m is met. Otherwise, at 
T > T m the FC forms, leading to a dramatic change 
of the Sommerfeld ratio 7(T). Indeed, according to the 
foregoing analysis, the major part of the FC momentum 
distribution is T-indcpcndent and therefore does not con- 
tribute to 7(T). Consequently, when formation of the FC 
is completed, the enhancement factor in 7(T) disappears 
in spite the huge FC density of states. 



B. Entropy excess 

A drastic change in the behavior of the entropy S(T) 
occur o 1 ' 16 ' 23 when, at T > T m , a FC is formed in the 
domain C under rearrangement of the momentum dis- 
tribution n(p) into the self-consistent FC solution n*(p). 
The basic entropy formula ([3]) of the original quasiparti- 
cle formalism remains intact: 



- / / [n*(p) lnn*(p)+(l— n*(p)) ln(l— n*(p))]du. 

(78) 
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However, due to the NFL component in n*(p), the system 
is seen to possess a T -independent entropy excess (p) . 
The situation we now face — with the strongly correlated 
fermion system having a finite value S* of the entropy 
at T — > resembles that encountered in a system of 
localized spins. In the spin system, the entropy referred 
to one spin is simply In 2, while in the system having a 
FC, we have S*/N ~ ?7ln2, where r\ — (pf — Pi)/pF is 
the dimensionless FC parameter. 

On the other hand, numerical calculations demonstrate 
that within the FC domain, the momentum distribution 
n*(p) changes rapidly under variation of the total den- 
sity p. The corresponding nonzero value of the derivative 
dS*/dp produces a huge enhancement of the thermal ex- 
pansion coefficient (3 ~ dS(T)/dp ~ n with respect to 
its FL value, proportional to T. 78 Consonant with this 
result, in many heavy-fermion metals it is found that 
/3 is indeed almost temperature-independent and that it 
exceeds typical values for ordinary metals by a factor 

1 3_ 10 4^9 



Furthermore, in the whole temperature interval from 
T = to T > T m , the spin susceptibility of a Fermi sys- 
tem beyond the QCP possesses Curie- Weiss-like behavior 
x{T) ~ 1/(T — Tw) with a negative Weiss temperature 
Tw- Measurements in 3 He films on various substrates 
and in numerous heavy-fermion compounds provide ex- 
amples of this NFL behavior. We emphasize that in our 
scenario, the negative sign of T\y holds even if the spin- 
spin interaction is attractive. This contrasts with the 
behavior in a system of localized spins, where the Weiss 
temperature has a negative sign only if the spin-spin in- 
teraction is repulsive. Moreover, in the case of localized 
spins, the Stoner factor must be suppressed, whereas in 
the vicinity of the QCP, this factor is usually enhanced. 

Another conspicuous feature of the physics beyond 
the QCP is associated with the Sommerfeld- Wilson ratio 
R sw = T X {T)/p 2 e C(T) ~ S*/C{T). Since the excess en- 
tropy 5* does not depend on T, it makes no contribution 
to the specific heat C(T); consequently in systems with 
a FC, one should see a great enhancement of i?sw- 



C. Curie- Weiss behavior of the spin susceptibility 



IX. CLASSICAL BEHAVIOR NEAR THE QCP 



Another peculiar feature of strongly correlated Fermi 
systems in the QCP region involves the temperature de- 
pendence of the spin susceptibility x{T) = Xo{T)/(l + 
<7 n (T)), where 



Xo = P 2 MT) = -2p\ 



2§ J n(p,T)(l 



dn(p, T) 
de(p) 



dv 



n(p, T))dv 



(79) 



and go is the spin-spin component of the interaction func- 
tion. As mentioned before, in the QCP region, the stan- 
dard FL structure is preserved only at T < T m . In par- 
ticular, the standard Pauli behavior x(T) = const, is 
maintained until "melting" of the multi-sheet structure 
occurs at T ~ T m , giving rise to flattening of the single- 
particle spectrum associated with a FC. 

At T > T m , insertion of n.*(p) into Eq. ([75]) yields the 
Curie-like term 



X*(T)=p 



2 Ceff(p) 



T 



with an effective Curie constant 



(80) 



(81) 



that depends dramatically on the density 23 ! 80 Since C e g 
is proportional to the FC parameter rj, we infer that 



Ces(p) = 2 / n»(p)(l - n*(p))dv 



a 



eff 



(82) 



Thus, all compounds in which the spin susceptibility ex- 
hibits the Curie-like behavior possess a large entropy. 



In previous sections, we have discussed the structure 
of the extended quasiparticle picture near the QCP and 
beyond it. In this section, we highlight one of the most 
distinctive (and counterintuitive) features of strongly cor- 
related Fermi systems in this domain of the Lifshitz phase 
diagram, notably their classical behavior. This aspect 
was first revealed in measurements of the specific heat 
C(T) of dense 3 Hc films at the lowest temperatures T ~ 1 
niK reached experimentally. As the temperature is low- 
ered into this regime, the specific specific heat C(T) be- 
haves a a 11 ' 12 



C(T)=/3 + 7 T, 



(83) 



thus exhibiting an unexpected T- independent term /3(p). 
Such behavior contrasts sharply with that of its 3D coun- 
terpart, bulk liquid 3 He, for which C(T) obeys standard 
FL theory in approaching zero linearly in T, with a pref- 
actor 7 proportional to the density of states N(0) oc M*. 
Indeed, FL theory works better and better as T is lowered 
toward zero, provided no supcrfluid transition intervenes. 

It is has been commonly assumec—i that the term /3(p) 
is related to the nature of the substrate which supports 
the 2D 3 He film. More specifically, it is thought that 
due to weak heterogeneity of the substrate (steps and 
edges on its surface), quasiparticles are delocalized from 
it, giving rise to the j3 feature of the heat capacity. 12 
This explanation is undermined somewhat by the fact 
that if substrate-induced disorder were responsible for 
the (3 feature, one would expect the departure of C(T) 
from the FL prediction to decrease as the film density in- 
creases, since the effects of disorder are most pronounced 
in weakly-interacting systems. Contrariwise, the specific- 
heat anomaly makes its appearance in just that den- 
sity region where the effective mass M* is enhanced and 
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2D liquid 3 He system becomes strongly correlated ! 11 ' 12 
In this situation, the impact of disorder should be sup- 
pressed. The weight of these considerations compells us 
to treat the observed behavior of C(T) as an intrinsic 
property of 2D liquid 3 He. 

Curiously, due to the presence of the residual term j3, 
this behavior mimics the classical Dulong-Petit (DP) Law 
of the specific heats C(T) of solids, the major contribu- 
tion to which comes from phonons, being obtained from 
the formula 



C B {T) oc / u(k) 



^dnsik) kdk 
~df 2tT' 



(84) 



where ns(fc) = [e^ k ^/ T — 1] _1 is the Bose-Einstein 
phonon distribution function. The T-dependence of 
Cb{T) is determined by relation between the temper- 
ature T and the Debye temperature Qd — u(fc max ) ~ 
cfcmax- At T > flu, the specific heat becomes indepen- 
dent of T, and the DP behavior is recovered. The greater 
the value of the parameter fc max characterizing the cut- 
off of the phonon spectrum, the more pronounced is the 
classical contribution (3 to the specific heat C(T), which 
is proportional to fc^ax for a system of dimensionality D. 

Normally, is sufficiently high that, true to its em- 
pirical discovery, the DP Law belongs to the realm of 
classical physics. However, we will see that in strongly 
correlated Fermi systems, the characteristic frequency 
of some collective mode can be extremely small such 
as to allow DP behavior (155)) in 2D liquid 3 He at mil- 
likelvin temperatures. Certainly, neither the usual hydro- 
dynamic sound nor longitudinal zero sound can qualify, 
since the group velocities of both modes remain finite as 
M* — > oo i 82 ' 83 The spin fluctuation mode is excluded as 
well, however, due to substantial Landau damping. The 
desired mode is in fact provided by transverse zero sound 
mode (TZSM), 84 whose properties near the QCP can be 
explicated in terms of its well-known dispersion relation 
in 3D matter^ 2 . 
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conveniently rewritten as 
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(86) 



with s = c/v F and Fx = fip F M*/ir 2 . The TZSM is seen 
to propagate only if F\ > 6, i.e., provided M* > 3M. 
Near the QCP where M*(p) — > oo, this requirement is 
always met. In this case, one has v F /c — >• 0, and Eq. 
simplifies to 
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0, (88) 



an analogous formula being obtained for a 2D system. 

It should now be clear that toward the QCP, the ef- 
fective Debye temperature fit = w(fc max ) = c^fc max goes 
down to zero, independently of the value of the wave num- 
ber fc max corresponding to the saturation of the TZSM 
spectrum. Thus, the necessary condition f2 t < T for 
emergence of a regime of classical behavior is always 
met. However, another condition must also be satisfied if 
there is to exist a well-pronounced classical domain at ex- 
tremely low temperature: the parameter fc max must not 
to be too small. 

There is no such a restriction in the conventional 
electron-phonon problem. Indeed, the phonon group ve- 
locity c depends weakly on the wave number k, whose 
characteristic value fc max coincides with its maximum 
possible value ~ l/^c the inverse distance between par- 
ticles in the system. As a result, the phonon contribution 
to the specific heat C(T) turns out to be proportional to 
the particle density p. However, softening of the TZSM is 
terminated around some rather small critical wave num- 
ber fc max , so that at greater wave numbers, u)(k) becomes 
larger than T, destroying classical behavior. Evaluation 
of the corresponding threshold value is a difficult prob- 
lem, because its solution requires analysis beyond the 
long-wave approximation. Here we consider the situation 
in systems with only a small proportion of FC, where this 
problem can be handled with some facility based on the 
familiar FL kinetic equatio n 39 ' 83 



(to — kv) 5n(p) = — kn 



dn{p) 
dp 



J r (p,Pi)5n(pi)dz;i. 



(89) 

Focusing on transverse zero sound in 2D liquid 3 He, we 
need only to include the term in the Landau amplitude 
J- involving the first harmonic f±. Making the usual 
identification (c t — cos 6)Sn(p) = (dn(p)/dp) <p(n), where 
cos 8 = kv/feu, Eq. (|89|) becomes 



= -flPFCOSt 



cos\ 



dn(pi)/dpi 
Ct — v(T) cost 



dpiddi 



" (27T) 2 ' 

(90) 

where cosx = cos^cos^i + sin 6 sin 9\, while v(T) is the 
FC group velocity, proportional to T. The solution de- 
scribing transverse zero sound is </>(n) ~ sin 9 cos 0. 

It is seen immediately that c t 3> v(T) ~ T; therefore 
the transverse sound in question does not suffer Landau 
damping. Then, upon retaining only the leading relevant 
term v(T) cos 9 /of of the expansion of l/(ct — v(T) cos 9) 
and performing straightforward manipulations, we are 
led to the simple result 
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(91) 



Factoring out an average value of the group velocity 
v(p,T) ~ T/pf, we find that 
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i.e., in the FC domain of the phase diagram, the velocity 
c t of the transverse mode depends on temperature so as 
to vanish like y/T as T — > 0. 

We now demonstrate that such a softening of the 
TZSM holds only as long as the wave number k does 



not exceed k n 



L 



Pf - Pi- 



The reason for this 



termination is that the noncondensed component of the 
quasiparticle system in 2D liquid 3 He comes into play 
at k > L, consisting of quasiparticles with normal T- 
independent dispersion de(p)/dp ~ pp/M. The group 
velocity c t (fc > L) of the transverse mode then soars up- 
ward, rendering the corresponding contribution to C(T) 
irrelevant. 

To evaluate this effect it is necessary to go beyond 
the long-wave approximation, a process beset with re- 
placement of kn9n(p) / dp by n(p + k) — n(p) and kv by 
e(p + k) — e(p). With these replacements, simple algebra 
leads to the behavior 

W 2 (fc) oc - J Hp) - n(p + k)) (e(p) - e(p + k)) 

(93) 

At k < L, almost all the FC states contribute to this 
expression on an equal footing, yielding relation (1921) . 
However, at k > L, the predominant contributions to 
the integral come from momentum regions where the 
difference |e(p) — e(p + k)| has its maximum value, which 
is T- independent. Thus we infer that at such large k val- 
ues, softening of the spectrum u>t(k) is terminated — in 
contrast to our previous claim 8 ^ that softening of ujt(k) 
persists until k > \fp~p~L. 

Since, as we have seen, k m: 
the result 



dv. 



L. Hence we arrive at 



^maxQ OC L 
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As long as the inequality L < (MT) 1//2 holds (or equiva- 
lently, T/e F > L/pp), the ratio Qt/T remains small, and 
the Dulong-Petit law C(T) = const, obtains. Thus, in 
spite of the low temperature, the specific heat behaves as 
if the system were situated in the classical regime. This 
paradoxical outcome is a consequence of the presence, in 
the strongly correlated Fermi system, of a macroscopic 
subsystem with heavy quasiparticles. As the tempera- 
ture ultimately goes down to zero at a fixed density p, 
the inequality L < (MT) 1 / 2 eventually fails, the quan- 
tum regime is restored, and the dominant contribution 
to the specific heat comes from the "normal" fermions. 
In other words, by reducing T sufficiently, one can reach 
a domain in which the FL behavior of the specific heat 
C(T) is recovered. 

Interestingly, the value of the constant term in C(T) 
can be evaluated in closed form in terms of the FC range 
L. Upon inserting Lu t (k) — ctk into Eq. (|84[) and inte- 
grating, the T-independent term in the specific heat is 
found to be 
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where N is the number of atoms in the film. 

We turn finally to a discussion of the impact of the 
TZSM on transport properties. In the case of a conven- 
tion Fermi liquid, the Fermi surface consists of a single 
sheet, so the TZSM has a single branch with velocity Ct 
exceeding the Fermi velocity vp. Consequently, emission 
and absorption of sound quanta by electrons is prohib- 
ited, and the role of the TZSM in kinetics is of little inter- 
est. However, in heavy-fermion metals, it is common for 
several bands to cross the Fermi surface simultaneously, 
thereby generating several zero-sound branches. For all 
branches but one the sound velocities are less than the 
largest Fermi velocity. Hence the aforementioned ban 
is lifted, and these branches of the TZSM spectrum do 
experience damping, in a situation similar to that for 
zero-spin sound. In the latter instance, Landau damping 
is so strong that the mode cannot propagate through the 
liquid. 82,83 It will be seen that this is not the case for 
damping of the TZSM, because of the softening of this 
mode close to the QCP. Due to the softening effect, the 
contribution of the damped TZSM to the collision inte- 
gral has the same form as the electron-phonon interaction 
at room temperature. 

To facilitate analysis of damping of the TZSM in sys- 
tems having a multi- connected Fermi surface, we restrict 
consideration to the case of two electron bands. The TPT 
is assumed to occur at one of the bands, so that its Fermi 
velocity, denoted again by Dp, tends to zero, while the 
Fermi velocity v of the other band remains unchanged 
through the critical density region. The model disper- 
sion relation for the complex sound velocity c = cr + ici 
becomes 
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It can easily be verified that the contribution of the 
second term to the real part of the right-hand side of 
Eq. is small compared to that of the first term, since 
vf/v q —> toward the QCP. On the other hand, not- 
ing that In [(cr + ici + v )/(cr + ici — v D )] ~ — in, the 
corresponding contribution mF^VpCn/ (4i> 2 ) to the imag- 
inary part of the right-hand side cannot be ignored, else 
cj = 0. By this reasoning, Eq. assumes the simpli- 
fied form 
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analogous to Eq. ([87]). Its solution obeys 



cr oc 



M M 
Cl oc 



M*(p) 



M*(p) 



(97) 



(98) 



Importantly, we see then that the ratio ci/cr oc 
yj M /M*{p) is suppressed in the QCP regime, which al- 
lows us to analyze the contribution of the TZSM to the 
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collision term entering the resistivity along the same lines 
as in the familiar case of the electron-phonon interaction. 



X. CONCLUSION 

Proceeding from the original FL quasiparticle picture 
due to Landau and Migdal, we have addressed the for- 
mation of flat bands in strongly correlated Fermi systems 
beyond a point where the necessary condition for stability 
of the Landau state is violated and hence subject to re- 
arrangement. Responding to a soaring exhortation from 
Migdal, ll Beri shireP (Embrace everything you can!), 
expressed to an interlocutor when A. B. approved what 
he heard, we have analyzed this phenomenon in diverse 
strongly correlated Fermi systems from neutron stars to 
atomic liquids to electron systems of solids. Absent, how- 
ever, are reviews of the latest achievements in the investi- 
gation of flattening of single-particle spectra in equations 
of particle physics^ in topological media ) 67 i 85 ' 86 includ- 
ing the analysis of flat bands on the surface of multi- 
layered graphene (see Refs. l87ll88l and references therein). 

We have seen that in the class of systems under con- 
sideration, rearrangement of the Landau state occurs by 
means of a cascade of topological transitions, at which 
the number of sheets of the Fermi surface grows steadily, 
as the spectrum of single-particle excitations e(p) ac- 
quires additional zeroes, in effect becoming flatter and 
flatter. Tracing the evolution of this spectrum as the 
coupling constant increases, we have shown that at zero 
temperature, the salient feature of the final stage of the 
evolution is the formation of the flat bands whose spec- 
trum e(p) proves to be dispersionless at T = 0, while at 
low T its dispersion becomes proportional to T. 

As we have seen in Sec. Ill, the underlying reason for 
such a universal rearrangement is based on the Pitaevskii 
identity; 40 i 41 derived from the Galilean invariance of the 
Hamiltonian of the system and gauge invariance, an iden- 
tity which coincides with the Landau equation for the 
spectrum e(p) as deduced from the assumption that e(p) 
is a functional of the quasiparticle momentum distribu- 
tion n(p). This coincidence establishes that the ground- 
state energy E is indeed a functional of n, because, ac- 
cording to the Lehmann expansion, e(p, n) is a variational 
derivative of E with respect to n. If one does not care 
about obedience to the Pauli principle, the global min- 
imum of the functional E(n) is then attained at some 
continuous function n*(p), devoid (as a rule) of jumps in 
momentum space. But we know that in weakly or moder- 
ately correlated Fermi systems, such solutions violate the 
Pauli restriction n(p) < 1, and, as per standard FL the- 
ory, the true quasiparticle momentum distribution turns 
out to be the step function 9(pf —p)- 

In strongly correlated Fermi systems, as we have seen 
in Sec. V.A, the situation changes drastically: begin- 
ning with a critical coupling constant Afc ; the inequality 
< n(p) < 1 is met, lifting the theoretical ban on the 
emergence of such an exotic creation as the fermion con- 



densate. 

Establishing this freedom does not in itself qualify as 
a rigorous proof of the viability of the swelling scenario, 
as such a proof must surmount obstacles associated with 
the interplay between single-particle and collective de- 
grees of freedom. The flattening of the single-particle 
spectrum may lead to softening of some spin/density- 
fluctuation mode, in principle generating a correspond- 
ing second-order phase transition before fermion conden- 
sation sets in. Alas, the detailed circumstances of this in- 
terplay await proper clarification. Thus it seems that we 
are forced to seek guidance from analysis of the available 
experimental data, even though, according to Migdal's 
criterion, these data are generally of secondary impor- 
tance. 

The flattening phenomena observed experimentally are 
so pervasive, however, that an exception to A.B.'s pol- 
icy statement is in order. The discovery of the quantum 
critical point, made almost simultaneously in 2D liquid 
3 He, MOSFETs, and heavy-fermion metals^^i 4 ! has 
been a milestone in experimental exploration of strongly 
correlated Fermi systems, providing a vigorous impe- 
tus to their theoretical investigation. Pivotal experi- 
mental guidance to theoretical development has emerged 
from evidence for the separation of QCPs from points of 
putative second-order phase transitions, uncovered first 
in 2D liquid He*2 and more recently in heavy-fermion 
metals. 56 Thus separation proves that collective degrees 
of freedom are, after all, not of crucial importance to the 
QCP phenomena under study. 

At this juncture in the development of QCP physics, 
one is prompted to ask whether we are facing a situa- 
tion similar to that surrounding the discovery of the W 
boson if in fact the phenomenon of flattening of single- 
particle spectra has already been observed in experiments 
on condensed-matter systems. If so, then the misinter- 
pretation of these experiments has impeded the revela- 
tion of a fundamentally new class of Fermi liquids. Un- 
fortunately, the scope of the data available for sufficiently 
definitive analysis is still quite limited. Measurements of 
angle-resolved photoemission electron spectra (ARPES) 
in solids, while otherwise promising, are not yet accu- 
rate enough to confirm or refute the linearity in T of the 
low-temperature dispersion of these spectra in relevant 
cases. 

Valuable insights might also be drawn from measure- 
ments of magnetic oscillations. In these measurements, 
the electron effective mass is extracted by performing a 
Lifshitz-Kosevich fit to the temperature dependence of 
the magnitude of observed oscillations. In strongly cor- 
related electron systems of high-T c superconductors and 
heavy-fermion metals possessing a QCP, such analysis 
is as yet limited to isolated examples. To our knowl- 
edge, a QCP has so far been documented only in the 
single high-T c superconductor YbBa2Cu306+a;, at a crit- 
ical doping x c ~ 0.52s where the measurements have 
been carried out at a temperature around 1 K in strong 
magnetic fields H > 55 T. Additionally, measurements^ 
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performed at much lower temperatures around 20 mK on 
the heavy- fermion metal CeCoIn 5 , which also has a QCP, 
have revealed deviations from the Lifshitz-Kosevich for- 
mula itself. Unfortunately, interpretation of the data on 
magnetic oscillations is burdened due to the necessity of 
imposing strong magnetic fields, which causes substan- 
tial distortion of the original electron motion, triggering 
magnetic breakdown^ 

In closing, let us briefly revisit the theoretical challenge 
presented by the growing body of thermodynamic mea- 
surements on strongly correlated Fermi systems. We fo- 
cus on the following striking features of the inferred ther- 
modynamic behavior at temperatures exceeding what 
may be interpreted as the critical temperature T m for 
melting of the FL structure and formation of a flat seg- 
ment of the quasiparticle spectrum e(p): 

(i) The existence of a T-independent entropy excess 
S*, reflected for example in a huge enhancement 
upon the FL value of the thermal expansion coeffi- 
cient /3 ~ /?* oc dS*/dP. 

(ii) A sharp falloff of the specific heat C(T) upon pass- 
ing through the critical temperature T m . 

(iii) Curie behavior of the magnetic susceptibility 
x(T) = C c ff/T, with a effective Curie constant dif- 
fering from the conventional value. 

We now comment specifically on each of these behav- 
iors in turn, within the context of the extension of the 
Landau-Migdal quasiparticle theory to embrace topolog- 
ical phase transitions, and especially fermion condensa- 
tion. 

(i) In ordinary FL theory, the entropy S is given by 
Eq. the curve S(T) starts at the origin and 
rises linearly with T. To attain entropy values 
S ~ In 2 at low T, the density of states, which 
specifies the slope dS(T)/dT, must be further en- 
hanced beyond the standard heavy-fermion boost 
associated with the width of the narrow /-band 
lying exactly at the Fermi surface. Nevertheless, 
experimental measurements 92 indicate that the en- 
tropy value S = 0.5 In 2 is realized in the heavy- 
fermion metal YbIr 2 SI 2 at T > T m ~ 1 K . More- 
over, another apparent manifestation of the flatten- 
ing phenomena is exhibited by the heavy-fermion 
metal CeCoIns mentioned above. Its thermal ex- 
pansion coefficient /3(T), measured in external mag- 
netic fields H ~ 5T imposed to suppress super- 
conductivity, increases strongly until T reaches the 
value T m = 0.3 Kp. whereas T > T m , /3(T) be- 
comes T-independent, being enhanced by a factor 
10 3 — 10 4 compared with standard FL values. 

(ii) A drop of the Sommerfeld ratio j(T) is observed 
in many strongly correlated electron systems as a 
temperature T m is exceeded. The record slump is 
found in YbIr2Si2 where the FL contribution to 



7(T) collapses at T > T m = 0.7 K. Usually such a 
slump is attributed to some second-order (antifer- 
romagnetic) phase transition. (As is well known, 
in the Landau theory of second-order phase tran- 
sitions there is a jump of the specific heat at the 
transition point.) Significantly, the experimental 
data on this compound and any other compounds 
having the QCP fail to show any jump of C(T) 
at any T. The scaling theory of second-order phase 
transitions also fails to explain the experimental be- 
havior of C(T). In addition, numerous attempts to 
establish the structure of the corresponding order 
parameter on the side of T m allegedly associated 
with the ordered phase, have also been unsuccess- 
ful. In light of all that has been done before in this 
article, none of this comes as a big surprise, and 
reminds us of a saying by Confucius: 

It is hard to find a black cat in a dark room, espe- 
cially if it is not there. 

Indeed, the arguments given in Sec. V.E explain the 
behavior of C(T) in the vicinity of T m , attribut- 
ing it to a crossover from a state with a multi- 
connected Fermi surface to a state having a flat 
band. Bearing in mind this association, we infer 
that in the QCP region, there are no hidden order 
parameters at all, since the behavior stems from 
topological phase transitions rather than second- 
order, symmetry-breaking phase transitions. 

(iii) Curie-like behavior of the magnetic susceptibility 
x(T > T rn ) = C c ff/T, which is incompatible with 
the customary FL Pauli behavior x(T) = const, 
was first observed in 2D liquid 3 He. 10 An analogous 
Curie-like behavior of x(T) has also been seen re- 
cently at very low temperatures in normal states of 
high-T c superconductors placed into a strong mag- 
netic field to suppress the superconductivit y 94 ! 95 
In both the cases, the effective Curie constant C e ff, 
being nontrivially dependent on the density p or 
doping x, is at variance with the ordinary Curie 
constant proportional to p. 

A detailed discussion of experimental support for or 
against the topological scenario for the quantum criti- 
cal point is beyond the scope of this article, especially 
bearing in mind that the measurements involved are for 
the most part very fresh, and — as cautioned by A. B. 
Migdal — may be subject to error or incomplete analysis. 
Thus, in summarizing the current state of knowledge, we 
are forced to recognize that the envisioned swelling of 
the Fermi surface in most strongly correlated condensed- 
matter many-fermion systems featuring the occurrence of 
the flat bands, has been neither validated nor disproved 
experimentally. The fermion condensate remains as elu- 
sive as the Cheshire Cat of Alice in Wonderland, teasing 
us with mischievous grins that become more and more 
visible and numerous. 
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